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ABSTRACT 

We present ultraviolet interstellar absorption line measurements for the sightline towards the 
09.5 V star /i Columbae (/ = 237?3, b = -27?1; d w 400 pc, z « 180 pc; (nm) w 0.06 cm^^) 
obtained with the Goddard High Resolution Spectrograph (GHRS) on board the Hubble Space 
Telescope. These archival data represent the most complete GHRS interstellar absorption line 
measurements for any line of sight towards an early-type star. The 3.5 km s~^ resolution of the 
instrument allow us to accurately derive the gas-phase column densities of many important ionic 
species in the diffuse warm neutral medium, including accounting for saturation effects in the 
data and for contamination from ionized gas along this sightline. We find that the effects of an 
H II region around n Col itself do not significantly affect our derivation of gas-phase abundances. 
For the low velocity material (—20 < wlsr ^ +15 km s~^) we use the apparent column density 
method to derive column densities. For the individual absorbing components at wlsr ~ —28.8, 
-1-20.1, -f31.0, and -1-41.2 km s~^, we apply component fitting techniques to derive column 
densities and 6-values. We have also used observations of interstellar Lya absorption taken 
with the GHRS intermediate resolution gratings to accurately derive the H I column density 
along this sightline. The resulting interstellar column density log A^(H I) = 19.86 ± 0.015 is in 
agreement with other determinations but is significantly more precise. 

The low-velocity material shows gas-phase abundance patterns similar to the warm cloud 
(cloud A) towards the disk star ^ Ophiuchi , while the component at wlsr ~ -1-20.1 km s~^ shows 
gas-phase abundances similar to those found in warm halo clouds. We find the velocity-integrated 
gas-phase abundances of Zn, P, and S relative to H along this sightline are indistinguishable 
from solar system abundances. We discuss the implications of our gas-phase abundance 
measurements for the composition of interstellar dust grains. We find a dust-phase abundance 
[(Fe -f Mg)/Si]^ = 2.7 — 3.3 in the low- velocity gas; therefore the dust cannot be composed solely 
of common silicate grains, but must also include oxides or pure iron grains. The low velocity 
material along this sightline is characterized by T « 6,000 — 7,000 K with ~ 0.3 cm"'^, 
derived from the ionization equilibrium of Mg and Ca. 

The relative ionic column density ratios of the intermediate velocity components at 
'^LSR = +31.0 and +41.2 km s^^ show the imprint both of elemental incorporation into 
grains and (photo)ionization. These clouds have low total hydrogen column densities 
(log A^(H) ~ 17.4 — 17.7), and our component fitting 6-values constrain the temperature in the 
highest velocity component to be T = 4, 000 ± 700 K. The electron density of this cloud is 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated bythe Association of Universities for Research in Astronomy, Inc. under the 
NASA contract NAS 5-26555. 
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rie ~ 0.6 cm~^, derived from the to '^P^j^ fine structure excitation of C II. The components 

at z^LSR ~ ~30 and —48 km s~^ along this sightHne hkely trace shocked gas with very low 
hydrogen column densities. The i^lsr ~ ~30 km s~^ component is detected in a few strong 
low-ionization lines, while both are easily detected in Si III. The relative column densities of the 
—30 km s~^ suggest the gas is collisionally ionized at moderate temperatures (T « 25,000 K). 
This is consistent with the measured 6- values of this component, though non-thermal motions 
likely contribute significantly to the observed breadths. 

Subject headings: dust, extinction - ISM: abundances - ISM: clouds - stars: individual 
(/U Columbae) - ultraviolet 

1. INTRODUCTION 

The measurement and analysis of interstellar absorption lines provides fundamental information on the 
content and physical conditions of the Galactic interstellar medium (ISM). In particular, measurements of 

the gas-phase abundances of the ISM have allowed us to infer the composition of interstellar dust and trace 
the variations in the dust make-up in a wide range of environments (e.g., Savage, Cardelli, & Sofia 1992; 
Spitzer & Fitzpatrick 1993, 1995; Sembach & Savage 1996; Lu et al. 1998). Furthermore, atomic absorption 
lines allow us to study the chemical evolutionary history of the Universe over 90% of its age through the 
study of gaseous QSO absorption lines (e.g., Lu et al. 1996; Pcttini et al. 1997, 1999; Prochaska & Wolfe 
1999), provided we can understand the imprint dust leaves on the measurements. The effects of dust are 
best constrained by our studies of gas-phase abundances at zero redshift. 

The study of Galactic interstellar absorption lines has been greatly aided by the Goddard High 
Resolution Spectrograph (GHRS) on board the Hubble Space Telescope (HST). The echelle-mode resolution 

of this instrument (FWHM « 3.5 km s^^) coupled with its ability to achieve high signal to noise ultraviolet 
observations of Galactic early-type stars makes it a very powerful instrument for studying absorption lines 
in the Galactic ISM. 

The star ^, Columbae is the GHRS high- and intermediate-resolution radiometric standard. As such, 
it was observed extensively during the Servicing Mission Orbital Verification (SMOV) stages after the 
installation of the Corrective Optics Space Telescope Axial Replacement (COSTAR). In the post-COSTAR 
era, more than 500 echelle-mode observations were made of /x Col using the GHRS, as well as a similar 
number of observations with the lower-resolution first-order gratings. This makes Col the most extensively 
observed early-type star with the GHRS. 

The design of the GHRS calibration observations during the SMOV period were such that the 

extensive /i Col cchcUc-modc datasct is characterized by extensive wavelength coverage and relatively high 
signal-to-noise ratios. The resulting high-resolution absorption line dataset is almost ideal for studying 
abundances along the low-density sightline to this star. Though this sightline has been extensively studied 
by the Copernicus satellite (Shull & York 1977; hereafter SY) and with the pre-COSTAR GHRS by Sofia, 
Savage & Cardelli (1993; hereafter SSC), the extremely rich dataset acquired after the installation of 
COSTAR represents a significant increase in resolution over the observations of SY and in S/N over those 
of SSC. The extensive wavelength coverage of the dataset includes observations of a wide assortment of 
ionic species, with most being observed in several transitions with a range of /-values. 

We have reduced and analyzed the extensive archival GHRS ultraviolet absorption line dataset for the 
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fi Col sightline. Our main objectives in this work are to very accurately derive the gas-phase elemental 
abundances in the low- and intermediate-velocity gas along this sightline as well as information on the 
physical conditions of the gas. From the gas-phase abundances we infer the dust content of the ISM in this 
direction and discuss the implications of these data for understanding the make-up of dust grains in the 
diffuse ISM. 

Our work is presented as follows. We discuss the properties of the fj. Col sightline and the previous 
studies of this sightline in §|[ In §|| we describe our reductions of the GHRS dataset, as well as our extraction 
and analysis methods for analyzing the ISM absorption line data. In §^ we discuss the gas-phase abundances 
of the low-velocity material along this sightline, including a detailed analysis of the contributions from an 
H II region about fi Col itself and the derivation of physical conditions for the absorbing gas. We present 
an analysis of the abundances of the intermediate velocity gas along this sightline in A discussion of 
the implications of this interstellar absorption line dataset is included in and we summarize our major 
conclusions in 



2. THE n COL SIGHTLINE 

The star fi Columbae (HD 38666) lies in the direction {I, b) = (237°, —27°) along a relatively 
unreddened sighthne [F,{B-V) = 0.02; Bastiaansen 1992]. Classified as an 09.5 V star by Walborn (1973), 
/i Col is a runaway star with a radial velocity relative to the local standard of rest (LSR) of wlsr = +90 
km s^^ (Gies 1987; Keenan & Dufton 1983). Lesh (1968) has placed the star as far away as 1000 pc. The 
absolute magnitude scale of Vacca, Garmany, & ShuU (1996) yields a slightly smaller distance of w 790 
pc, given the observed magnitude V = 5.16. However, recent Hipparcos measurements have constrained 
its parallax to be tt = 2.52 ± 0.55 milliarcseconds (Ferryman et al. 1997), implying a distance of only 

pc. This result is more consistent with the distance derived from Stromgren photometry (410 pc) 
by Keenan & Dufton (1983). This discrepancy in distance scale may be related to problems in assigning 
luminosity classes to rapidly rotating spectral standards (Lamers et al. 1997). We will adopt the Hipparcos 
results throughout this work, though the reader should be aware of the continuing gaps in our knowledge of 
the stellar distance scale. 

The adopted distance, coupled with the observed neutral hydrogen column density of 
logA^(Hl) = 19.86 ±0.015 (see Appendix g), implies an average line of sight density of (nn) ~ 0.06 
cm~'^. Molecular hydrogen makes a negligible contribution to the total hydrogen column density with 
log A^(H2) = 15.50 summed over J = to 4 (Spitzer, Cochran, & Hirshfeld 1975). 

Highly ionized gas along the ^ Col sightline implies the presence of both hot (T > 10^ K) collisionally 
ionized material and photoionized gas near the star. York (1974) fit the strong interstellar O VI absorption 
along this sightline with a column density log A^(0 VI) = 13.82 ± 0.02 and a Doppler parameter 6ovi = 43.8 
km s~^. Brandt et al. (1999) present GHRS observations of weak C IV, Si IV and N V absorption at 
3.5 km s'^^ resolution. Figures show the normalized interstellar absorption line profiles for the /i Col 
sightline observed with the GHRS, including the profiles of C IV and Si IV (see Figure^). The profile widths 
of the highly ionized atoms increase with ionization potential. The breadth of the profile implies T < 2 x 10^ 
K, while the width of the C IV line, for example, implies temperatures T < 2.3 x 10''' K; much of the velocity 
widths may be due to non-thermal motions. The O VI and C IV profiles are centered near (wlsr) = ^2 
km s"\ while the Si IV is centered at (wlsr) = +5 km s"^ The ratio iV(C IV)/iV(0 VI) = 0.11 ± 0.01 
is consistent with that observed for other sightlines intercepting Galactic disk gas (Spitzer 1996). Brandt 
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et al. suggest much of the O VI and some C IV arises in an evolved supernova remnant along the line of 
sight, though some of the high- ion absorption likely arises in the interface between the Local Cloud and the 
Local Bubble. They also show that a significant amount of the Si IV column density is likely produced in 
a low-density H II region surrounding /i Col. We will discuss in more detail the possible contributions of 
ionized gas to the absorption line measurements. 

The low-ionization material along the sightline to /i Col was studied with the Copernicus satellite by 
SY at a resolution of « 13 km s~^ and with the echelle gratings of the pre-COSTAR GHRS by SSC with 
a resolution of « 3.5 km s~^. These studies have identified low-ion absorption in four main absorbing 
components, the properties of which are given in Table |^. This table includes the number by which we 
will refer to each of these absorbing regions, the approximate central velocity and the range of velocity 
over which the absorption from each region extends, both in the LSR frame, as well as the identifications 
of these regions in the nomenclature of SY and SSC. Regions 2, 3, and 4 are relatively well separated in 
velocity and may represent individual absorbing "clouds;" region 1 is a blend of several clouds that overlap 
in velocity. The region we designate as component 5, which appears at negative velocities, was described by 
SY as "trailing absorption." Each of the absorbing regions in Table |l| shows different gas-phase abundances, 
suggesting changing patterns of elemental incorporation into dust grains and ionization. In general the 
absorption along this sightline is consistent with the Routly-Spitzer effect (Routly & Spitzer 1952), where 
the gas-phase abundances of refractory elements increase with velocity (SY; SSC; Hobbs 1978). Also given 
in Table |^ are the temperatures, T, and electron densities, rie, for each component from our analysis below 
(see and sfj). 

Lockman (1991) has discussed the H I 21-cm emission profile towards Col, which is reproduced in 
Figure ||. These data, from Lockman, Hobbs, & ShuU (1986), have a 21' beam and have been corrected for 
stray radiation. The column of H I from the 21-cm data is more than twice that derived from observations 
of Lyman-a along the sightline to /i Col, implying much of the H I column observed in 21-cm radiation 
comes from beyond the star. In the direction of ji Col Galactic rotation is expected to carry the gas to 
positive velocities, with a maximum velocity of vlsb. ~ +5.2 ± 0.6 km s^^. The majority of the H I emission 
and the UV-absorption (i.e., the region 1 blend) resides at velocities allowed by Galactic rotation. A quite 
narrow component at wlsr ~ 3 km s~^ shows the peak brightness temperature in the H I emission profile, 
though several broad absorbing clouds also contribute to the blend of region 1. It is clear that the H I 
content of components 2, 3, and 4 is significantly less than that of component 1. Though UV absorption 
profiles along this sightline show gas out to wlsr ~ +42 km s~^, the H I emission profile shows an extended 
wing to wlsr ^ +60 km s"^ 

Component 1 is centered at wlsr ~ +3 km s^^ and contains ~ 90% of the neutral gas along the 
sightline, as evidenced by the relative strengths of the absorbing regions in the undepleted species S II. 
Ground based observations of the Na I D2 line at a resolution of 1 km s^^ show the complex to consist 
of several (at least three or four) velocity components (Hobbs 1978; see Figure ||). Though the absolute 
velocity of the absorption is uncertain, SY tentatively associate the observed H2 absorption with component 
1. 

The papers of SY and SSC have shown absorbing complex 1 likely has depletion characteristics 
similar to that of the warm cloud towards C, Oph (component A of Savage et al. 1992). The low velocity 
resolution of the Copernicus data make the separation of absorption from components 1 and 2 difficult; the 
higher resolution observations of SSC are superior in this respect. Those elements with high condensation 
temperatures, such as Fe II and Cr II, exhibit gas phase abundances that are sub-solar with respect to S II 
by w 1.3 dex. Sofia et al. (SSC) suggest that the level of depletion exhibited by the /i Col component 1 and 
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the warm component along the C Oph sighthne may be typical of the low-density, warm neutral medium 
(WNM) of the Galactic disk. 

The presence of ionized material in the velocity range encompassed by component 1 is suggested by 
N II and Si III absorption in the data of SY. The fine structure lines of N II**, which trace the densest 
ionized regions, are centered at vi^sn = to +2 km s^^. ShuU & York argue that the majority of the ionized 
gas along the sightline is associated with an H II region surrounding fi Col and estimate a density (rig) « 0.2 
cm~'^. Recent observations of this sightline with the WHAM Fabry-Perot spectrometer (Reynolds et al. 
1998), which are reproduced in Figure |8[ show that the ionized gas is indeed concentrated near ulsr = 
(M. Haffner, 1998, priv. comm.). The imprint of ionized regions on the column densities of (primarily) 
neutral gas tracers can be one of the largest uncertainties in studying the gas-phase abundances in the 



Galactic WNM, such as in component 1. We will discuss this contamination of our dataset in §4.1 



Component 2 appears distinctly in species that tend to be highly depleted in component 1. A quick 
comparison of the Ca II and Na I profiles of Hobbs (1978) shows that the ratio of these two species changes 
significantly as the velocity of the gas increases (these data are reproduced in Figure ||). This difference 
is likely due to the return of elements to the gas phase due to dust destruction. ShuU & York (SY) and 
Sofia et al. (SSC) have shown that the refractory elements have a much higher relative abundance in this 
component than in the lower velocity material. This cloud is at a velocity that would place it well beyond 
the star if it were simply participating in Galactic rotation. 

The intermediate-velocity gas along the /i Col sightline has been less-well observed due to the relatively 
low column of material present in these clouds. Centered at ulsr ~ +33 and +42 km s~^, respectively, 
components 3 and 4 show patterns of abundances quite different than the low-velocity material. Both 
SY and SSC examine the gas-phase abundances in component 4, though component 3 was only identified 
by SSC. In both cases, the effects of ionization may be significant, though the gas-phase abundances of 
these clouds may also be high, suggesting substantial grain destruction (SSC). ShuU & York also noted the 
presence of "trailing" absorption in the Si III profile. This material, which we identify as component 5, has 
a velocity relative to the local standard of rest vlsr ~ —30 km s^^, which is inconsistent with Galactic 
rotation. 

Some 90% of the gas towards /i Col is associated with component 1, which has velocities roughly 
consistent with expectations due to Galactic rotation. Most of the remaining 10% of the material has 
velocities which are inconsistent with rotation. Sembach & Danks (1994) have found that on average ~ 10% 
of Ca II absorption is at forbidden velocities. They estimate a cloud-to-cloud velocity dispersion in this 
forbidden- velocity gas of cr « 22 km s^^. The forbidden- velocity gas towards /i Col, therefore, seems not to 
show highly unusual kinematics compared with the observations of low-density sightlines by Sembach & 
Danks. 



3. DATA PROCESSING AND ANALYSIS 

Table ^ lists basic information of the individual spectra used in our analysis of the /i Col sightline. 
The present set of observations was acquired for the purpose of evaluating the in-flight performance 
of and flux-calibrating the high-resolution modes of the GHRS following the installation of COSTArJ^ 



^Details about the GHRS and its in-flight performance characteristics can be found in Robinson et al. (1998) and Heap et 
al. (1995). 
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The data have been collected over a large span of time, beginning in early 1994 after the installation of 
COSTAR. As such, the observations do not represent a completely homogeneous data set. We rely most 
heavily on measurements made with the echelle-mode Ech-A and Ech-B gratings, which give a resolution of 
« 3.5 km s^^ (FWHM). An extensive dataset exists for the first-order gratings as well, and we have used 
observations taken with the G140M and G160M gratings to derive the column densities of H I (see Appendix 
^) and Fe III along the fi Col sightline. Typically the star was observed for 30-120 seconds through the 
large science aperture (LSA; l'.'74xl'.'74) with four substeps per diode. Appropriate measurements were 
made of the inter-order scattered light in all cases (see Cardelli et al. 1993), and the observations employed 
the comb-addition routine with the on-board doppler compensator enabled. We have in general restricted 
ourselves to using the post-COSTAR data for this sightline due to the degradation in resolution of the 
pre-COSTAR LSA observations. We have, however, made use of pre-COSTAR small science aperture (SSA; 
0'.'22x0'.'22) data in our component fitting analysis (see § ^.3.2| ). 

3.1. Data Processing 

Our calibration and reduction of the data follows procedures similar to those discussed in Savage et al. 
(1992) and Cardelli et al. (1995). Our determination and propagation of errors follows Sembach & Savage 
(1992) for our measurements of the integrated equivalent widths and column densities. The basic calibration 
was performed at the GHRS computing facility at the Goddard Space Flight Center and at the University 
of Wisconsin-Madison using the standard CALHRS routine.^ The CALHRS processing includes conversion of 
raw counts to count rates and corrections for particle radiation contamination, dark counts, known diode 
nonuniformities, paired pulse events and scattered light. The wavelength calibration was derived from the 
standard calibration tables and should be accurate to approximately ±3.5 km s^^. 

The final data reduction was performed using software developed and tested at the University of 
Wisconsin-Madison. This includes the merging of individual spectra and allowing for additional refinements 
to the scattered light correction. The inter-order scattered light removal discussed by Cardelli et al. (1990, 

1993) is based upon extensive pre-flight and in-orbit analysis of GHRS data and is used by the CALHRS 
routine; the coefficients derived by these authors are appropriate for observations made through the SSA. 
We find that many of the LSA observations required an additional correction to bring the cores of strongly 
saturated lines to the appropriate zero level. The final scattered light coefficients, dc, used for each group 
of spectra are given in Table ^j. Many of the observations had no strong lines which would allow us to 
refine the values of dc] those cases for which we have adopted the Cardelli et al. (1993) values, having no 
independent measure of dc , are marked with a colon in Table |^. In general the signal-to-noise ratio and 
spacing of the individual spectra did not warrant solving explicitly for the FPN spectrum. Those regions 
for which we have derived the noise spectrum and removed it (following the algorithm of Cardelli & Ebbets 

1994) are identified in Table |. 

To bring all of the species into a common velocity reference we have applied a "bootstrap" technique 
similar to that discussed in Cardelli et al. (1995). We have aligned in velocity space lines of the same species 
found in different observations. We have further attempted to align several ions with similar ionization 
and depletion characteristics, or similar velocity structure. For example, the strong lines of Mg I and II, 
Fe II, Al II, and Si II have been shifted in velocity space to be brought into alignment with the Si II A1304 



^CALHRS is part of the standard Space Telescope Science Institute pipeline and the STSDAS IRAF reduction package. It is 
also distributed via the GHRS Instrument Definition Team for the IDL package. 
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line. Component 4 is well separated from lower-velocity absorption for the stronger lines of these species, 
making the alignment relatively straight-forward. The Si II A1808 line was brought into this reference 
frame using the wings of Si II A1304, and the weak lines of Fe II have been aligned with the wings and the 
distinct component 2 of the stronger Fe II lines. The O I A1302 line appears in the same exposures as Si II 
A1304 and was used to tie N I and C II into the velocity frame of the more depleted species. The lines of 
the heavily depleted species Ni II A1370 and Cr II AA2056, 2062 were aligned to the weaker AA2249 and 
2261 lines of Fe II; the other lines of these species lacked sufficient signal-to-noise to enable any correction 
to their velocity zero-points. Cr II A2062 was then used to align Zn II A2062, which appears in the same 
observations. The latter line was then used to bring Zn II A2026, P II A1152, and S II A1250 into the 
common velocity reference. The two stronger lines of S II were aligned to the A1250 line. Lastly, the lines 
Mg II AA1239 and 1240, and Mn II AA2577 and 2594 were aligned as well as possible with Si II A1808. The 
Mg II to Si II alignment should be reasonably secure given the somewhat similar component structure of 
these transitions, but the Mn II lines show different component structure, particularly for component 2 
centered near ulsr ~ +21 km s~^, which seems to be indicative of the different depletion characteristics of 
these elements. The alignment of these last lines is more uncertain than most. 

To determine an absolute velocity frame we have measured the heliocentric velocity of the O I* A1304 
and O I** A1306 telluric absorption lines, which are present in the same specrum as the Si IIA1304 line. We 
have determined the central velocities of each of these telluric lines and compared those with the velocity 
of the spacecraft at the time of the observation plus the correction for the Earth's motion towards the 
star. A correction of —0.3 km s~^ was needed to bring the Si II A1304 line observed through the SSA into 
the heliocentric rest-frame. An additional correction was then made to convert heliocentric velocities to 
the LSR frame. Assuming a solar neighborhood speed of -1-16.5 km s^^ in the direction {l,b) = (53°, 25°) 
(Mihalas & Binney 1981) implies ulsr — Vheiio = —16.5 km s^^ for the sightline to /i Col. We have, however, 
applied a shift of —19.9 km s^^ to all heliocentric velocities in order to be consistent with previous studies 
of this sightline which have assumed a solar motion of -1-20 km s^^ in the direction {l,b) = (56°, 22°) 
[{a, 5)i9oo = (18'', -^30°); see York & Rogerson 1976; also adopted by York 1974, SY, and SSC]. 



3.2. Absorption Profiles and Measurements 

Continuum normalized interstellar line profiles for all species treated in this work are shown in Figures 
|l|— 0. Each profile was normalized by fitting low-order (< 5) Legendre polynomials to the local stellar 
continuum in regions free from interstellar absorption (Sembach & Savage 1992). In general the continuum 
of the star, which has a radial velocity wlsr — +90 km s~^ (Keenan & Dufton 1983) and projected rotational 
velocity vsini w 111 km s^^ (Penny 1996), was well behaved, making the fit to the stellar continuum 
relatively certain. In some cases, however, the interstellar absorption coincides with stellar lines in a way 
that makes the continuum placement more ambiguous. Examples of such occurences include the lines Si II 
A1193, Si III A1206, and N I A1199. We have marked lines with less than certain continuum placement in 
the tables of data presented herein. For comparison with the low-ionization GHRS data presented here, we 
also include in Figure H] ground based absorption profiles of Ti II from Welsh et al. (1997) and the Ca II and 
Na I profiles from Hobbs (1978) as well as the H I emission profile from Lockman, Hobbs, & ShuU (1986). 
Also shown in Figure || is the WHAM spectrum of Ha emission along this sightline (M. Haffner, 1998, priv. 
comm.). The GHRS data for Si IV and C IV from Brandt et al. (1999) are shown in Figure 

The integrated equivalent widths, Wx, are given in Table ^ for each species, along with the la error 



estimates (see §3.3.1). The range over which the equivalent width and apparent column density integrations 
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extend for each absorbing region are given in Table Also listed in Table || are the ionization potentials of 
the measured ionic species and the next lower ionization state of the element, vacuum wavelengths, adopted 
values of the oscillator strengths for each transition and the studies from which we have drawn these values, 
and the empirically estimated signal-to-noise ratios. Oscillator strengths used in our analysis are generally 
taken from the compilation of Morton (1991), using the recommended updates listed in Table 2 of Savage 
& Sembach (1996a) with a few updates for new determinations detailed below. 

For the Ni II and the weak Mg II transitions recent determinations suggest oscillator strength revisions 
by factors of ~ 2. We choose to adopt the recent empirical determination of the oscillator strengths for 
Mg II AA1239 and 1240 by Fitzpatrick (1997). These /-values are determined from a comparison of the 
strong Mg II transitions near A2800 with those near A1240 in high signal to noise GHRS observations. 
Fitzpatrick's recommended /-values are a factor of ~ 2.4 larger than the Hibbert et al. (1983) theoretical 
calculation and a factor of ~ 2 smaller than the emperical determination of Sofia, Cardelli, & Savage (1994). 

The Ni II oscillator strengths are derived from a combination of the Fedchak & Lawler (1999) and 
Zsargo & Federman (1998) results. Zsargo & Federman (1998) have placed many of the Ni II /-values 
on a consistent relative scale using GHRS observations of several stars. Their compilation includes all of 
the transitions we observe with the exception of the A1317 line. We find no significant evidence in our 
data that the ratios of the AA1317 and 1370 /-values from Morton (1991) should be modified, though the 
signal-to-noise ratio for the latter line is less than ideal for this type of study. Fadchek & Lawler (1999) 
have very recently provided absolute laboratory measurements of the oscillator strengths of a number of 
vacuum ultraviolet Ni II transitions, including the AA1709 and 1741 transitions observed in this work. The 
ratios of the absolute /-values derived by these authors for the AA1709 and 1741 transitions are in excellent 
agreement with the results of Zsargo & Federman. Fedchak & Lawler suggest using the /-values derived 
by Zsargo & Federman multiplied by a scale factor of 0.534 ± 0.05. We adopt this recommendation in this 
work, using values log A/ that are —0.272 dex below the values suggested by Zsargo & Federman. We note 
however that there is a discrepancy between the implied column densities derived from the A1741 transition 
and the AA1317 and 1370 transitions. The measurements of the A1741 transition in our dataset are based 
on a single Ech-B exposure with S/N ~ 21. The possibility exists that a FPN feature is present in these 
data; however, the velocity structure appears similar to the weak Cr II and Fe II lines with comparable 
signal-to-noise ratios. We believe that there may still be uncertainties in the relative oscillator strengths 
between the lines longward and shortward of 1700 A in our dataset. The change in oscillator strengths 
suggested by Fedchak & Lawler (1999) and Zsargo & Federman (1998) not only has implications for the 
dust content of diffuse interstellar clouds in the Milky Way but also, and perhaps more importantly, clouds 
the interpretation of [Ni/Fe] measurements in high-redshift damped Lya systems (e.g., Lu et al. 1996; 
Prochaska & Wolfe 1999; Kulkarni, Fall, & Truran 1997). We discuss the implications of the new Ni II 
oscillator strengths in more detail in §|^. 

Our estimation of the errors inherent to our measurements of the integrated equivalent widths and 
apparent column densities (see below) includes contributions from photon statistics, continuum placement 
uncertainties, and zero- level uncertainties (Appendix A of Sembach & Savage 1992). We have adopted a 
2% zero-level uncertainty throughout. Though this may overestimate the errors in regions near heavily 
saturated lines, we feel it is appropriate given the uncertain scattered-light properties of the LSA. 
Continuum placement uncertainties were estimated based upon the effects of adjusting the continuum level 
by ±0.4 times the rms noise about the fit. These sources of error are independent and have been added in 
quadrature to produce the final error estimate quoted with our measurements. 

The sources of error discussed above make no allowance for the existence of FPN features in our data. 



- 9 - 



The strength of FPN features is reduced significantly by co-adding spectra that have been shifted along the 
diode array from one another, so that features constant in diode-space are shifted in wavelength-space. Even 
given this improvement, FPN features are found in our reduced data, often mimicking weak interstellar 
absorption lines. For example a FPN feature is present at wlsr ~ +41 km s^^ in our data for the S II 
A1253 line (see Figure p, almost exactly coincident with the expected absorption from component 4. In 
this case, however, we can identify it as a FPN feature because it appears in only one of the two co-added 
observations of this wavelength region. For those wavelength regions covered by only one observation, our 
ability to discriminate between true interstellar absorption and weak FPN features becomes less robust. In 
these cases we are aided by the excellent wavelength coverage of the current /i Col dataset: we are often 
able to check the reality of absorption features in many different transitions of the same ionic species. 

The nominal short-wavelength limit of the GHRS is 1150 A, given the inefhciency of the magnesium 
fluoride coatings of the _ff5'T optics at wavelengths less than 1150 A. However, the short wavelength Digicon 
detector has a LiF window, and observations shortward of this are possible. We include in our data 
measurements of the Fe III A1122.526 transition as well as the N I triplet at A1134, and transitions of Fe II 
at AA1133, 1143, and 1145. Unfortunately, two of the ground-state transitions of neutral carbon are nearly 
coincident with the Fe III transition, at AA1122.518 and 1122.438, which lie at —2.1 and —23.5 km s~^ 
relative to the Fe III velocity zero point, respectively. We do not believe the C I contamination of Fe III is 
a significant problem for the following reason. We detect the much stronger C I line at A1560.309 with an 
integrated equivalent width of Wx = 14.7 ± 0.9 mA. This is equivalent to a combined equivalent width from 
the C I transitions at AA1122.518 and 1122.438 of Wx = 0.66 ± 0.04 mA, implying the C I transitions make 
a negligible contribution to the Fe III measurement {Wx = 19 ± 4 mA). 

3.3. Analysis Methods 

In theory one can separate an observed interstellar absorption profile into individual absorbing clouds 
along the line of sight, if the data fully resolve these individual entities. By fitting models for the absorption 
from each cloud, or component, one can determine the column densities, central velocities, and Doppler 
parameters for each absorbing cloud along a given line of sight. An excellent example of applying this 
approach to study the abundances in the diffuse ISM is the work of Spitzer & Fitzpatrick (1993), who used 
the GHRS to study the abundances and physical conditions towards HD 93521. In practice one can run into 
significant uncertainties with this approach, particularly for clouds closely-spaced in velocity. An example 
of a region where the component fitting techniques become difficult, leading to a lack of uniqueness, is 
the principal absorbing region along the line of sight towards /i Col (component 1; —17.0 ^ VhSR 15.5 
km s~^). The Fe II and Mg I profiles suggest that there are multiple blended components in this velocity 
range. However, constraining a fit to the data for the less depleted species of S II and Si II is more difficult. 

To derive accurate column densities for the components along the line of sight to /i Col, we will apply 
the component fitting technique for the higher velocity gas (components 2-5) but will primarily rely upon 
the apparent column density, Na{v), method described by Savage & Sembach (1991) for dealing with the 
central low- velocity blend we have designated component 1 . 
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3.3.1. Apparent Column Density Method 

In analyzing the low-velocity components towards Col, we will make use of the so-called apparent 
column density, or Na{v), method, which gives information on the velocity structure of the absorbing 
material that is model-independent (Savage & Sembach 1991). In short, a continuum normalized absorption 
profile I{v) = e~'^'^^'"\ for a transition having wavelength. A, and an oscillator strength, /, is related to the 
apparent column density per unit velocity, Na{v), by 

A^» = ^^ = 3.768 x 10-^, (1) 
Tre^ /A /A(A) 

in miits atoms cm^^ (km s^^)~^, where A is given in A. In the absence of unresolved saturated structure, 
which can be identified by comparing 7Va(w)-profiles for different transitions of the same species. Savage & 
Sembach (1991) have shown this method provides a valid, instrumentally blurred representation of the true 
column density as a function of velocity, N{v). 

Examples of Na{v) profiles for the (presumably) non-depleted species S II, the moderately-depleted 
Mg II and Mn II, and the highly-depleted Fe II, are shown in Figure |[ For each of these ionic species two 
transitions with different oscillator strengths are plotted. One can see that the examples we have chosen, 
in general, show good agreement between the two transitions. An example where our profiles exhibit 
unresolved saturated structure is seen in the Fe II profiles. The Na{v) values near -yLSR ~ +21 km s^^ 
are lower in the stronger A2586 line than the weaker A2374 line. This is evidence for unresolved saturated 
structure in component 2 within these Fe II lines. For Fe II, however, it is still possible to accurately derive 
N{v) with transitions weaker than the A2586 transition. 

Table ^ contains the velocity-integrated apparent column densities (Savage & Sembach 1991) and 
estimated errors for the transitions studied in this work. The sources of errors for these integrations were 
taken to be the same as those described above for deriving the integrated equivalent widths. We have not 
included the uncertainties in the /-values in this error budget. In the absence of unresolved saturated 
structure, these column densities are equivalent to the true column densities in the velocity ranges outlined 
in Table |l| The column densities are representative of the individual components considered here if there 
exists no significant blending between the components. 



3.3.2. Component Fitting 

The column densities derived through the Na{v) method above may be subject to large uncertainties 
in cases where the individual absorbing regions overlap significantly in velocity. In particular component 2 
may be heavily contaminated by overlap with absorption due to component 1 in lightly depleted species 
(cf.. Figure 0). To more cleanly separate the higher velocity gas seen in components 2, 3, 4, and 5 we use 
component fitting techniques to derive the interstellar column densities. 

The component fitting approach begins with a model of the interstellar absorption spectrum, which 
consists of k individual components described by their central velocities, Vk, Doppler spread parameters, 
bk, and column densities, N^. The individual model components are assumed to be well approximated by a 
Voigt profile. This model is then convolved with an appropriate instrumental line spread function (LSF), 
and the value of minimized between this blurred model and the observed line profile to determine the 
best-fit parameters. Our component fitting analysis makes use of software kindly provided by E. Fitzpatrick 
(1998, priv. comm.) and described in Spitzer & Fitzpatrick (1993) and Fitzpatrick & Spitzer (1997). 
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Where they exist, we have used observations taken through the SSA of the GHRS in our component 
fitting analysis. While these data often have lower signal to noise, the LSF of the SSA has been carefully 
studied by Spitzer & Fitzpatrick (1993) and has a slightly better resolution than the LSF of the LSA. 
However, most of the data presented here were taken through the GHRS LSA. The LSF for the GHRS LSA 
has not been well characterized. Robinson et al. (1998) present a LSF for the post-COSTAR LSA; however, 
we found this LSF inadequately matched the results derived from our fitting of SSA profiles and was not 
able to fit narrow deep lines correctly (e.g., component 4). In Appendix ^ we derive a new LSF for the 
post-COSTAR GHRS LSA. The new LSF is a sum of a strong narrow Gaussian and a weak broad Gaussian. 
The narrow component has a FWHM of 1.09 diodes, while the weak component has a FWHM of ~ 4.21 
diodes with a peak approximately 4.5% that of the narrow component at A2800. Therefore the broad weak 
component contains 15.1% of the spread function area. This fraction is a function of wavelength, however. 
We discuss this LSF in more detail in Appendix 

The best fit 6-values and column densities from our component fitting analysis are given in Table |. 
The best fit central velocities are given in Table We fit the blend making up component 1 with three 
components having approximate central velocities (wlsr) ~ —8, 1, and 7 km s~^, though there is some 
variation in the best fit values. Due to the lack of uniqueness in this central blend, we do not report here 
the results for these individual components, but only the sum of their column densities in Table |^. Our 
purpose in fitting component 1 was not to disentangle this blend, but to approximately account for the 
overlap of this region with the more distinct component 2. 

Where several lines exist for a given ionic species, we have fit all of the profiles simultaneously. For 
several species it was necessary to adopt 6- values derived from fits to other ions. This approximation was 
necessary when either the signal-to-noise ratio of the spectra were not high enough for the fitting to be 
reliable (e.g., Cr II or Ni II) or when the profiles of the components were not distinct enough to provide 
appropriate information to constrain the fit (e.g,. for component 2 in P II and Zn II). In these cases we 
have adopted relatively well-constrained ^-values from lines of similar atomic mass. To assess the error 
contribution of the adopted &- values to our derived column densities, we have also calculated the best fit 
models for h ± la. The differences between the bzLla results and the best fit 6-value results were added in 
quadrature to the formal fitting error. For the S II profile we have fixed the central velocity of component 2 
to that derived for Si II. This was done because the unconstrained fit yielded velocity structure and relative 
column densities in the central blend that were significantly different than that of any other ion, and these 
differences impacted the fitted component 2 parameters. We found, however, that by holding the velocity 
of this component we were able to obtain a fit in good agreement with our results for the other ions. 

In general the component fitting results for component 1 agree with the Na{v) integrations, suggesting 
little in the way of unresolved saturated structure or confusion from component overlap. We find the 
column densities derived for component 2 in lightly depleted species are typically ~ 0.2 dex lower than the 
results derived from integrating the Na{v) profiles. This is a result of the overlap from wings of components 
in the low- velocity blend that are included in the Na{v) integration. We see that component 4 may be 
significantly saturated in a number of profiles by comparing the component fitting results with the column 
densities derived from a straight integration of the Na{v) profiles (e.g., Mg II and O I). The &-values derived 
from our component fitting analysis are b < 2.4 km s~^ (FWHM^ 4.0 km s^^). Thus it is important to use 
the component fitting results for this cloud. 
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3.4. Adopted Column Densities 

To derive the best column densities for component 1, we take the weighted average of all the transitions 
for a given species that show no evidence for unresolved saturated structure in their Na{v) profiles. We 
present our adopted final column densities in Table ^ noting where we have chosen not to use an observed 
transition in deriving these column densities. We are relatively certain that the column densities in Table 
for the low-velocity absorbing components toward /i Col are not significantly affected by saturation effects. 
The individual clouds making up the blend of component 1 are relatively broad, and the absorption profile 
seems to be fully resolved by the echelle-mode resolution of 3.5 km s^^. 

In a few cases we have chosen to use the component fitting results for component 1. These cases have 
been marked in Table and include Ni II, Fe II, and Zn II. The Zn II observations show a systematic offset 
of 0.08 dex between the integrated column densities of the two transitions. This cannot be due to saturation 
effects since the stronger of the two lines gives a higher apparent column density, exactly the opposite of 
the expected behavior in the presence of saturation. This behavior may be caused by uncertainties in the 
oscillator strengths of the transitions. We have chosen to adopt the component fitting results for Zn II 
because we believe the column densities are a better compromise between the two profiles than the Na{v) 
values and the formal errors are more representative of the true errors. 

The final adopted column densities for components 2 — 5 are from our component fitting results. Most 
of the values given in Table are the result of simultaneously fitting all the available transitions of a given 
species, with exceptions noted in the table. Table |^ also includes our derivation of the H I column density 
along this sightline, which is described in Appendix 0. 



4. ABUNDANCES IN THE LOW VELOCITY GAS 

In this section we discuss the observed abundances in the low-velocity absorbing regions (components 
1 and 2) along the line of sight to /i Col and the implications of these abundances for interstellar dust. This 
velocity range (from wlsr ~ to +29 km s^^) not only contains the vast majority of the warm neutral 
absorbing column but also material associated with ionized gas along the path to fi Col. It is necessary 
to examine the degree to which material in primarily ionized gas may affect the derivation of relative 
abundances in the neutral material along this sightline. We discuss this contamination and our assessment 



of it in §4.1 



Throughout this paper we will be discussing the normalized gas-phase abundance of elemental species. 
We define the normalized gas-phase elemental abundance of a species X relative to F as a function of 
velocity to be 

[X/Y], EE \og{Na{v)x/Na{v)Y} - log{X/Y}Q, (2) 

where Na{v)x and Na{v)Y are the apparent column density per unit velocity for the elements X and 
Y, respectively, and the quantity {X/Y}q is the solar or cosmic reference abundance ratio of the two 
species X and Y (e.g., the meteoritic abundances from Anders & Grevesse 1989). We will use [^Z^] to 
denote the equivalent of Equation (H) when one uses velocity-integrated total column densities in place of 
Na{v)x and Na{v)Y- This nomenclature is equivalent to others' definition of the logarithmic depletion, 
D{X) (e.g., SSC; Spitzer & Fitzpatrick 1993). When deriving gas-phase abundances it is typically assumed 
w [X/y], if X+ and are the dominant stages of ionization in the warm neutral medium. For 
the most part this assumption is justified, though the effects of ionized gas along the line of sight may 
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modify the ions in a different manner, causing this assumption to break down. We will show in §4.1 that 
this is not a significant effect for component 1. However, for components 3 — 5 this assumption will not be 
appropriate (see §^). 



4.1. Ionization Effects 

The presence of absorption due to the ions S III, Si III, Si IV, Al III, and Fe III in our GHRS spectra 
and the Copernicus observations of N II strongly suggest the presence of ionized hydrogen (H"*') along the 
line of sight to n Col. The WHAM observations of this region imply the presence of ionized gas in this 
direction at velocities compatable with those of the ionized tracers observed by GHRS, though much of the 
emission may come from beyond the star. For regions primarily containing H°, an element X whose first 
ionization potential falls below that of hydrogen is predominantly found in its singly ionized form, . 
Thus measurements of X+/H° or X'^ /Y^ are generally good indicators of the gas-phase abundance of 
the element X in the neutral material. However, the inferred presence of H+ along the sightline to /i Col 
complicates this simple picture since the relative contributions of X+ and may be different for each 
element in the H+-containing region and will be dependent upon the ionization structure of the region. 
Therefore, it is important to investigate the effects of H+-containing regions along the line of sight on our 
derived gas-phase abundances. 

/i Col is an 09.5 V star with T^fj ~ 33, 000 K (Howarth & Prinja 1989)0 and is therefore hot enough 
to ionize its immediate surroundings. A SIMBAD search of the area within 5° of /i Col reveals two B2.5 IV 
stars along the sighthne to /x Col: 7 Col and HD 41534. These stars, with Tefj ~ 20,000 K, he ^20 pc 
from the sightline to the star and should contribute very little to the ionized column along the sightline 
compared with /i Col itself. Within 10° of the star there are a total of seven B stars of types B2.5 or later. 
Given the lack of any 0-type stars within 300 pc of the /z Col sightline (SY) and the late spectral types of 
the B-type stars found near the sightline, we will continue under the assumption that the majority of the 
warm ionized gas that occurs along this sightline is contained in a photoionized nebula around /x Col. The 
column of H II along the line of sight can be predicted using 

7V(H+) = iV(H° + H+) -iV(H°) (3) 
« iV(S+) • {H/S}0 - iV(H°). 

Using this approach SY predict log A^(H+) = 19.94. Using our accurate determination of A^(H I) towards 
fi Col (Appendix^, coupled with our high-quality S+ measurements, we predict logA^(H+) = 19.17 ±0.14. 
This estimate relies on the key assumption that the total sulfur to hydrogen abundance along this sightline 
is solar. The WHAM spectra and the GHRS observations of S III, Al III and Si 11* (see below) suggest that 
ionized gas contamination is most significant for component 1. 

To more reliably assess the contribution of material associated with an ionized nebula about /i Col to 
our absorption line measurements, we have used the photoionization equilibrium code CLOUDY (v90.04; 
Ferland 1996, Ferland et al. 1998) to model the ionization and temperature structure of such an H II 
region. To estimate the stellar spectrum from /i Col for use in our photoionization models, we adopt 



^There are varying determinations of the stellar effective temperature. Keenan & Dufton (1983) have published T^ff ~ 
31,400 K, while the temperature scale of Vacca et al. (1996) suggests Tf,ff 34,600 K is appropriate. We will adopt the 
intermediate Howarth & Prinja (1989) value. Changes over this range of temperatures do not significantly affect the results of 
this section. 
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an ATLAS line-blanketed model atmosphere (Kurucz 1991) with the stellar parameters Tg// — 33,000, 
\og{g) = 4.0 and L* = 2.6 x 10^ Lq = 1 x 10'^^ ergs s"^, close to the values estimated by Howarth & 
Prinja (1989) for ^l Coif] For a range of densities and filling factors of the ambient ISM, we have calculated 
the ionization and temperature structure of model nebulae from a distance 0.03 pc from the star to the 
point where the electron density, Ue, falls to 5% of the total hydrogen density, rin. We have used solar 
abundances throughout]^ The ambient densities used in the models presented here are nn = 0.02,0.05, 0.2, 
and 0.5 cm~^. The densities 0.05 and 0.2 cm"'^ are approximately the average line of sight density, and 



the estimated density of the H II region gas from the excited states of N II (SY) and Si II (sec §4.4). The 
highest value was used to study regions more dense than the limit of Ue ^ 0.2 cm~^, while the lowest is 
used to show the relative constancy of the results. The first three models are those presented in Brandt et 
al. (1999) to discuss the source of the Si IV absorption along the line of sight to fx Col. Similar models have 
also been presented by Howk & Savage (1999) to derive the gas-phase abundance of [Al/S] in the ionized 
medium of the Galaxy. 

In order to match the observational constraints, we assume all of the column density of S III arises in 
the photoionized region about the star. For the densities considered, we vary the volume filling factor of 
the material until a match to the observed S III column density of logiV(S III) = 13.82 is obtained. Table ^ 
contains the physical parameters for each of our models and the predicted column densities of important 
ionic species. The predicted column densities of H+ for the models given in Table ^ are in the range 
log A^(H II) = 19.01 — 19.10, which is in rough agreement with the H II column density predicted by scaling 
our S II column density according to Equation (^). The H II column densities from our CLOUDY models 
are mildly sensitive to where the model nebulae are truncated. 

Among the most important conclusions we have drawn from our modelling of the ^ Col H II region is 
that the ratio of Af(Al+^)/7V(S+2) predicted to occur in an H II region is relatively insensitive to the model 
assumptions. Howk & Savage (1999) have used similar models and observed values of A^(A1^^) and A^(S+^) 
to determine the normalized gas-phase abundance of Al to S in the H II region surrounding /it Col. Their 
model-corrected result, 

[A1/S]hii = -0.78 ±0.08, 

is insensitive to the adopted density or filling factor of the nebula, and the ionizing photon flux of the star. 
Further, the sensitivity to the shape of the ionizing spectrum (i.e., to T^ff), is also relatively small. The 
error estimate comes from the standard deviation about the mean predicted value of iV(Al+^)/iV(S+2) for 
models using input spectra in the range 27,000 < Te// < 39,000 K (see Howk & Savage 1999). 

A similar calculation can be made regarding the gas-phase abundance of Fe relative to S in the H II 
region about /x Col. Howk & Savage (1999) have shown that values of 7V(Fe'''^)/iV(S"'"^) predicted by the 
models are much more sensitive to uncertainties in the input effective temperature of the ionizing source. 
Using observations of A''(Fe~'"^) and N{S^^) and CLOUDY models, they determine 

[Fe/S]Hii = -0.87 ±0.21 



®As noted above, some of the fundamental properties of this star have a range of pubUshed values. We continue to adopt 
the intermediate Howarth & Prinja (1989) results, but see also the results of Keenan & Dufton (1983) and the relationships 
with spectral types given in Vacca et al. (1996). 

^We have run models with Orion nebula abundances and abundances appropriate for warm disk gas to assess the effects of 
the different abundances on the temperature and ionization structure of the nebula. We have also included opacity due to dust 
grains to test the robustness of our results. The derived results from these models are completely consistent with our approach. 
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in the /i Col H II region. 

The sub-solar abundances of Al and Fe in the /i Col H II region likely implies the existence of dust in 
the ionized ISM about this star. An Al depletion of 0.8 dex is similar to that found for refractory elements 
along warm disk+halo sightlines (Sembach & Savage 1996). The similar depletion of Fe is also consistent 
with warm disk+halo sightlines. 

For the undepleted ions S II and P II, which are the dominant ionization stages in the WNM, our 
models predict absorption from material in the H II region may account for ^ 10% of the total observed 
absorbing column density. Thus the presence of a low-density H II region surrounding fi Col may add a 
systematic contribution of « 0.04 to 0.05 dex to the column densities of undepleted species in component 1. 
This is a relatively small contribution. We have chosen not to report the photoionization model results for 
the undepleted species Zn II given the large uncertainties in the adopted atomic parameters, particularly 
the recombination coefRcients and ionization cross sections, for this element. 

Assessing the impact of the H II region to the measurements of highly depleted species is more 
complicated. We have assumed solar abundances for our model H II region. Given our analysis of the Al III 
and Fe III absorption above, using solar abundance models will over-estimate the contribution of the H II 
region to the total line of sight absorption measurements for depleted species. This can easily be seen in 
the results tabulated for Fe II in Table |[ In some cases we predict more Fe II absorption than is actually 
measured for component 1. Because there is evidence for sub-solar abundances in the measured ratios 
iV(Fe+^)/iV(S+2) and N{Al+'^)/N{S+^), we have confidence that dust exists in the ionized gas about fi Col. 

The estimated abundance [Fe/S] in the H II region is roughly consistent with that of warm disk+halo 
sightlines. We can use this abundance to estimate the corresponding abundances of the elements Mg, Si, 
and Mn in the H II region. The compiled logarithmic depletions (equivalent to our normalized gas-phase 
abundances) of these elements are tabulated in Savage & Sembach (1996) for such sightlines to be 

[Mg/H] « -0.31; [Si/H] « -0.25; and [Mn/H] « -0.66, 

where we have corrected the value of [Mg/H] for our adopted /-values. If we apply these depletions to our 
model H II region results, the relative contribution of H II region gas to the measured total column densities 
of these species is approximately 

Mg II : 0.05 dex; Si II : 0.06 dex; Mn II : 0.09 dex; and Fe II : 0.06 dex, 

where we have assumed [Fe/H]= —0.87. Thus, material associated with the H II region may provide a 
relatively small, though not insignificant, contribution (uncertainty) to the measured column densities. This 
contribution is not large enough to hinder our major conclusions regarding the gas-phase abundances of the 
primarily neutral medium. We will assume that the contribution of H II region gas to the measurements of 
Zn II is similar to that of the undepleted elements S II and P II. For Cr II and Ni II we find results similar 
to those of Fe II and Mn II, respectively. It is important to point out that the expected contribution of 
an ionized nebula to the depleted and non-depleted species are very similar. Thus while the contributions 
to the S II and Fe II column densities from the ionized nebula are of order 0.05 to 0.06 dex, the ratio 
A^(Fe II)/iV(S II) is only changed by ~ 0.01 dex. However, when comparing singly-ionized species to neutral 
hydrogen, this contribution is more significant. 

The aforementioned uncertainties in the fundamental stellar parameters (temperature, luminosity, and 
surface gravity) do not seriously affect our estimates of contamination from gas in an H II region around 
fi Col. The temperature and surface gravity of the star change the shape of the spectrum, but over the 
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range of allowable values, we have found varying these parameters changes little in our models. Changing 
the luminosity does not affect the results at all. In this case the important parameter is the "ionization 
parameter" (see Howk & Savage 1999), which can be made made constant with varying luminosity by 
simply adjusting the ambient density and/or filling factor. 

It is important to constrain the ionization fraction of the primarily neutral gas along the sightline. 
There may be partially ionized gas in the H I-containing regions along the sightline that do not contribute 
to the S III used to constrain our H II region models, but which may contribute (particularly) to the S II 
column density along this sightline. We can roughly estimate the ionization fraction = ne/nn of the 
neutral regions using the Ar I measurements of SY with the work of Sofia & Jenkins (1998; hereafter SJ). 
Sofia & Jenkins have suggested that the intrinsic Ar/H abundance in the ISM is close to to the solar system 
abundance [following SJ we adopt 12.0 + log(Ar/H)0 = 6.52], and that sub-solar Ar I/H I measurements 
reflect the over-ionization of Ar I relative to H I. This over-ionization of Ar I is a result of its very large 
ionization cross section relative to that of H I. Sofia & Jenkins develop a formalism to relate the observed 
discrepancy, [Ar/H], between the observed and expected Ar abundance to the ionization fraction, Xe, along 
an interstellar sightline (see their Equation 11). Their treatment incorporates the relevant ionization and 
recombination rates, effects of charge exchange reactions, and photoionization into high ionization states of 
Ar through a quantity they denote P'ai- Characteristic values of P'ai for various interstellar conditions 
and ionizing spectra are tabulated in their Table 4. Using the Ar I measurements of SY, corrected for the 
suggested Ar I oscillator strengths used in SJ, we find [Ar/H] = —0.24 ± 0.10 towards Col using our 
iV(H I) measurement. Thus it would seem that some degree of ionization is present even in the neutral gas 
along this sightline, unless the reference abundance of Ar is incorrect. Using Equation (11) from SJ, with an 
estimate of P'ai = 10 ± 5 from their Table 4, we find Xe = 0.09 ± 0.07. Therefore we expect the contribution 
of ionized hydrogen to the total neutral column density due to partially-ionized gas to be 0.037 ± 0.028 dex. 
Thus there could be an additional correction of ^ —0.04 dex needed to account for partially-ionized gas in 
the H I regions along the /i Col sightline. 

We note that the relative abundance [Ar/S] derived with Ar I and S II measurements is extremely 
sensitive to relatively small amounts of ionized gas along a given sightline since the measured gas phase 
S II and Ar I increase and decrease, respectively, as the level of ionization increases. Furthermore, neither 
of these elements is expected to be incorporated into dust grains in large amounts (see SJ), and their 
nucleosynthetic origins are the same, both being a-process elements primarily produced in high-mass stars. 
The usefulness of [Ar/H] measurements for identifying the existence of partially-ionized gas has been 
emphasized by SJ; we suggest that the [Ar/S] ratio, which can be determined as a function of velocity, is 
an even more sensitive indicator of the combined effects of partially-ionized gas in H I-bearing regions and 
H II region contamination. 



4.2. Gas-Phase Abundances 

Table ^ contains the values of [X/S] in the low-velocity components 1 and 2 for the elements considered 
here. For comparison with other sightlines and QSO absorption line systems we also give the sightline 
integrated values of [^/H] for each element X. These results are derived from the adopted column 
densities presented in Table 0. Also given in this table is the adopted solar reference abundance for each 
element relative to hydrogen, and the normalized gas-phase abundances for components 3 and 4, which 
will be discussed in The quoted errors here contain only those described in §3.2. Our adopted solar 
reference abundance system is taken from Savage & Sembach (1996a), who rely primarily on the meteoritic 
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abundances of Anders & Grevesse (1989) with updates for C, N, and O by Grevesse & Noels (1993). 

We have plotted the data from Table ^ in Figure |l^. The ordinate of these plots is the gas-phase 
abundance for the elements listed along the abscissa, normalized to S. The top panel of Figure [l0| contains 
only the data from components 1 and 2. The bottom panel also shows these data, but we have overlayed 
the values [X/S] for the warm cloud along the ( Oph sightline (cloud A of Savage et al. 1992) as the dashed 
line, and for the spread of warm halo cloud abundances from Sembach & Savage (1996) as the hatched 
region. Data for the gas-phase abundances of interstellar Ti in individual halo clouds are sparse. For Ti 
we plot an upper limit equal to that of Ni and a lower limit that is equal to the lower limit of Fe in these 
clouds. The values of [X/S] from Savage et al. (1992) and Sembach & Savage (1996) have been adjusted to 
reflect our choice of oscillator strengths. 

The values of [AT/S] for elements in the component 1 blend are quite similar to those found in the 
warm ^ Oph cloud. This was also the conclusion of SSC, though with considerably more uncertainty. This 
level of sub-solar abundances is consistent with other "warm disk" -like clouds as categorized by Sembach & 
Savage (1996). Similarly, the relative abundance pattern seen for component 2 is mostly within the small 
range of "halo" cloud values. This is somewhat deceiving, however, since the earlier GHRS data for the 
/i Col sightline (SSC) were used by Sembach & Sembach (1996) in defining this spread of values. The new 
data presented here for component 2 support the claim of Sembach & Savage that the upper-limit of [A/S] 
for many elements shows relatively little cloud-to-cloud variation about the mean for each species. What is 
perhaps surprising, given the new distance to /z Col, is that a cloud within w 400 pc of the sun {z <^ 180 
pc) has abundance patterns similar to clouds at relatively large heights from the plane. If one adopts the 
distance derived from the Vacca et al. (1996) absolute magnitude scale, namely d ~ 790 pc, the z-height 
becomes z ~ 360 pc from the midplane. This is more than one H I scale height above the midplane. 

In the gas of component 1, the values [Si/S] and [Mg/S] are slightly higher than those seen along the 
C Oph sighthne warm cloud. Using the oscillator strengths suggested by Fitzpatrick (1997), Si and Mg 
show quite similar abundance patterns. Given the similarity of the Si II A1808 line profile to those of Mg II 
AA1239, 1240, this is not a surprising result. The similar abundance behavior of these elements is seen even 
more clearly in plots of [A/S]i,, the gas-phase abundance of an element A as a function of velocity. 



Figure 11 presents plots of the normalized gas-phase abundances for several elements as a function of 
velocity, using the adopted solar system abundances listed in Table ^ (principally the meteoritic abundances 
of Anders & Grevesse 1989). We have plotted the ratios of P II, Zn II, Mg II, Si II, Mn II, Fe II, Ni II, and 
Cr II to S II in this figure. For all of the species but Ni II, the width of the data bins is half a resolution 
element (i.e., ~ 1.75 km s^^ or two of the original data points for profiles taken with four substeps per 
diode). We have binned the Ni II profile to one point per resolution element. The error bars represent the 
sources of errors discussed in §^ as well as a contribution from possible velocity scale offsets. The latter 
was derived by shifting one profile relative to the other by half a resolution element (in both the positive and 



negative velocity directions) and adding the resulting errors in quadrature to those of §3.2. These velocity 
shift error estimates are the primary cause of the asymmetric error bars in Figure |l^. The range over which 
data points are plotted in this figure is determined by the significance of each point once we have added all 
of the sources of noise. This presentation assumes that the Na{v) w N{v) and [A+/S+]t, « [A/S]t,. The 
former of these assumptions is justified by our examination of the Na{v) profiles of each of the ions for 
which unresolved saturation might be present. In no case do we find evidence for saturation effects in the 
transitions presented in Figure ^ The latter of these assumptions is likely valid for material in the warm 
neutral medium. In the lower-right of each panel we have included a bar representing the maximum degree 
of uncertainty the /i Col H II region is thought to add to the individual measurements, as discussed in the 
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previous subsection. 

There is an increase in the gas-phase abundances of the elements Fe, Mn, Cr, Ni and possibly Si near 
''^LSR ~ +20 km s~^, which is associated with the significant strengthening of component 2 in moderately- 
and heavily-depleted species. However, one of the more striking aspects of Figure |l^ is the relative 
constancy of the ratios plotted here over the range —10 < wlsr ^ +14 km s~^. The absorbing clouds 
making up the component 1 blend have similar gas-phase abundances. The [X/S]„ profiles for Si and Mg 
are quite similar, with [Si/S]„ « [Mg/S]u « —0.4 dex. The Mg abundance profile is slightly fiatter than the 
Si profile. Figures |l| and|ll| show that Mg and Si trace each other well (see Fitzpatrick 1997). 

The abundances [X/S]u for the elements P and Zn are consistent with solar system abundances, 
i.e., [P/S]t, « [Zn/S]„ ~ 0. This behavior holds over the whole velocity range considered in Figure |l^. 
This supports our choice of solar system abundances as a reference for this dataset. B stars in the solar 
neighborhood may have lower intrinsic abundances than the sun (e.g., Kilian-Montenbruck et al. 1994; Gies 
& Lambert 1992), suggesting the proper "cosmic" reference abundance may be sub-solar by ^ 0.2 dex. In 
the case of the sightline towards fi Col, we observe gas with solar system abundance ratios of P and Zn to 
S, which we believe argues for using solar system reference abundances for these clouds. The choice of solar 
relative abundances for these clouds is particularly strong given the sources of metal production for S, an 
a-element, are thought to be different stars than those producing Zn, which traces Fe-peak elements. Given 
these elemental abundance ratios, we will rely primarily upon the solar system abundances of Anders & 
Grevesse (1989) as a cosmic reference (as compiled in Savage & Sembach 1996a). 



4.3. Implications for Interstellar Dust 

The sub-solar gas-phase abundance patterns seen in Figures ^ and ^ likely represent the imprint of 



elemental incorporation into interstellar dust grains. Our arguments from §4.1 suggest that the ionized gas 
contribution to the measured total column of each of the species considered in Figures |l^ and|ll| IS mmor, 
certainly less than ^0.1 dex. The depletion of elements from the gas- to the solid-phase is a well known 
phenomenon and allows us to infer the elemental make-up of dust grains in the diffuse ISM (Savage & 
Sembach 1996a). 

The striking similarities between the normalized abundances in the warm ( Oph cloud and the material 
making up component 1 along the ^ Col sightline suggest that this level of depletion, or incorporation of 
elements into the dust-phase, is relatively common among low H I column density clouds that make up 
the warm neutral medium in the solar neighborhood (SSC). The H I column densities of these two cloud 
complexes are quite similar [log A^(H I)^; oph A ~ 19.74; Savage et al. 1992]. 

Of interest for determining the types of grains present in the ISM is the dust-phase abundance of each 
element with respect to hydrogen. The dust-phase abundance of a species X relative to hydrogen, (X/H)^, 
is given by 

(X/H), = (X/H),-(X/H)g, (4) 

where the subscripts d, c, and g refer to the dust-phase, cosmic, and gas-phase abundances of X, 
respectively. Table |l^ gives the dust-phase abundances of a series of elements relative to H (given in parts 
per million H) and relative to Si for components 1 and 2. We also give the dust-phase abundances of Mg, Si, 
and Fe when adopting B-star abundances as determined by Gies & Lambert (1992) and Kilian-Montenbruck 
et al. (1994). Gies & Lambert do not derive the abundance of Mg in their work, so we have adopted 
the Mg abundance from Kilian-Montenbruck et al. in this case. The values of {X/}i)c are given for each 
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element in the different reference systems. The values (X/H)^ have been derived assuming S is present in 
its cosmic abundance, i.e., is not depleted into grains. The cosmic sulfur abundances in the three systems 
are: (S/H)^ = 1.9 x 10"^ (Anders & Grevesse 1989); (S/H)^ = 1.6 x 10"^ (Gies & Lambert 1992); and 
(S/H)c = 9.3 X 10^^ (Kihan-Montcnbruck et al. 1994). We have not included P or Zn in this table, as 
neither shows any evidence for incorporation into dust grains in our data. 

Among the elements considered here, the most abundant in dust grains are Mg, Si, and Fe. The solid 
forms of these elements in the ISM are thought to include silicates, oxides and possibly metallic iron. 
Among the silicate forms thought to be most common are various pyroxenes, (Mg, Fe)Si03, and olivines, 
(Mg, Fe)2Si04 (Ossenkopf et al. 1992). If the Mg- and Fe-bearing dust in the ISM towards /i Col were 
primarily made up of only these types of silicates, one would expect a ratio of [(Fe+Mg)/Si]£i « 1 — 2. 
Assuming the solar abundances given in Table fol we find: 



If one adopts the Kilian-Montcnbruck et al. (1994) B-star abundances as the cosmic reference, these values 
become: [(Fe + Mg)/Si]d = 5.06 ±0.24 in component 1; and [(Fe + Mg)/Si]<j = 15± 14 in component 2. The 
abundances of Gies & Lambert (1992) yield intermediate values. The composition of dust grains containing 
Mg, Si, and Fe in the material making up component 1 is inconsistent with the sole consituent of this dust 
being silicate pyroxenes or olivines, irrespective of the cosmic abundance one chooses to adopt. Component 
2 may be similar in this regard, but the large uncertainties make this statement much less certain. 

Whittet et al. (1997) have recently reported on the detection of O— H stretching modes in OH groups 
along the diffuse sighthne to Cygnus 0B2 No. 12 (VI Cygni No. 12) using the Infrared Space Observatory. 
Given the lack of ices along this sightline, they attribute this feature to hydrated silicates. The values of 
(Mg/Si)c( in components 1 and 2, assuming solar system abundances, is consistent with the incorporation of 
Mg into a mixture of the phyllosilicates talc, Mg3Si40io[OII]2, and serpentine, Mg3Si205[OII]4. However, 
Whittet et al. find a very small fraction of silicates along that sightline are hydrated. Thus we believe it 
unlikely phyllosilicates can provide for enough of the Mg to allow Fe-bearing silicates to account for the 
total dust-phase Fe abundance. Some amount of the Mg- and Fe-bearing dust is therefore likely in the 
form of oxides or pure iron grains. Examples of Mg-, and Fe-bearing oxides include MgO, FeO, Fe203, and 
Fe304 (Nuth & Hecht 1990; Fadeyev 1988). 

It is clear from Figures |l^ and |l^ that there is a larger gas-phase abundance of the refractory elements 
in component 2 than in component 1. In particular, the gas-phase abundances of Fe, Cr, and Ni show 
much higher abundances in Figure |l^. In principle higher gas-phase abundances of Fe-peak elements in 
interstellar clouds could be interpreted as evidence for enrichment by Type la SNe (Jenkins & Wallerstein 
1996). We discount this mechanism for providing the enhanced gas-phase abundances seen in component 
2 over component 1 along the /i Col sightline, suggesting instead that the large enhancements seen in 
the Fe-peak elements over the other elements in Figure |l^ are a result of the return to the gas phase of 
highly depleted elements. The increases in the gas-phase abundances per million H in component 2 over 
component 1 are: 15.4 ± 1.1 for Fe, 11 ± 4 for Si, 0.018 ± 0.006 for Ti, 0.27 ± 0.04 for Cr, 0.078 ± 0.015 for 
Mn, and 0.54 ± 0.07 for Ni. If we were to assume gas with abundances similar to component 1 had been 
enriched by gas from a Type la SN, we would expect factors of 2.5 and 3 — 7 more Mn and Ni, respectively, 
relative to Fe using the nucleosynthetic yields of Nomoto et al. (1984) and Thielemann et al. (1986). Also, 
relative to the models of Type la SN nucleosynthesis, the observed increase in Si relative to Fe is a factor of 
^ 2 too large, while the increase in Ti to Fe is too high by a factor of 2.5 to 160, depending upon the exact 



[(Fe + Mg)/Si], 
[(Fe + Mg)/Si], 



2.70 ±0.11 in component 1; and 
3.3 ± 1.4 in component 2. 
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nucleosynthesis result used. Therefore, we conclude that the increase in gas-phase abundances observed in 
component 2 relative to component 1 is likely not due to enrichment by the nucleosynthetic products of a 
Type la SN. 

We interpret the higher gas-phase abundances seen in component 2 as a return of elements to the gas 
phase from the solid phase in material that has been processed by a shock(s). The dust-phase abundances 
of many elements in component 2 have been lowered relative to those in component 1. Examining the data 
in Table ^ we see that relative to component 1, the dust- phase abundances of Si, Mg, and Fe have changed 
by approximately a factor of two. The values of [(Fe -I- Mg)/Si]d in component 2 are consistent with those of 
component 1 given the large errors for this component. There is evidence for lower dust-phase abundances 
of Si, Fe, Ti, Cr, Mn, and Ni in component 2. If the dust-phase abundances of the warm ( Oph cloud and 
the clouds making up component 1 are characteristic of a standard depletion in the low-density WNM of 
the Galaxy, we can estimate the fraction of material returned to the gas-phase from dust in component 2. 
The data in Table |l| suggest the return of (55 ± 20)% of Si, (50 ± 4)% of Fe, (23 ± 8)% of Ti, (60 ± 9)% of 
Cr, (26 ± 5)% of Mn, and (32 ± 4)% of Ni to the gas phase from dust grains, when compared to the values 
appropriate for component 1. 

If one interprets the differences in abundance patterns between components 1 and 2 as evidence for 
the stripping of grain mantles from standard dust-phase abundance of the WNM, the composition of the 
mantles may be inferred from a direct comparison of the dust-phase abundances of components 1 and 2. 
The data for the /x Col sightline then suggest that the mantles of dust grains are consistent with silicate 
and oxide components. Our large errors in the values for component 2 make a detailed determination of the 
mantle composition difficult. 



4.4. Physical Conditions 

In this subsection we discuss the information about the physical conditions in the low-velocity material 
towards fi Col contained in our absorption line data. The approaches used here are well-discussed in the 
works of Spitzer & Fitzpatrick (1993, 1995) and Fitzpatrick & Spitzer (1994, 1997). 



Electron Densities in the Ionized Gas 

Copernicus observations of N II, N 11*, and N II** have yielded estimates for the average electron density 
in the ionized gas, {ne,i), along the sightline to /i Col (SY). These authors' estimates yield (rig^i) ~ 0.16 
to 0.22 cm~'^ for the ionized gas. We detect weak absorption arising from the excited fine structure 

level of Si II at A1264.738 in our GHRS dataset. This allows us to verify the results of SY. The equation for 
collisional excitation equilibrium of the Si^ levels may be written 

v42l"-(Si+*) + 72l71e7l(Si+) = 7i2nen(Si+), (5) 

where 712 and 721 are the rate coefficient for collisional excitation and de-excitation, respectively, and 
A21 the spontaneous downward transition probability. In the case of ion excitation by electrons we write 
(Spitzer 1978) 

8.63 X 10-6 /-£;i2\ 3 -1 ff■^ 

T12 = g^ri/2 exp j cm s . (6) 
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For warm gas like that expected in an H II region, exp{Ei2 / kT) w 1.0 (to the accuracies being considered 
here). We adopt the values A21 = 2.17 x 10^^ s~^ (Mendoza 1983; his Appendix), fli2 = 5.58 (Osterbrock 
1989; his Table 3.3), and the statistical weight gi = 2. Replacing the particle densities with column densities 
we have 

, , i,,iV(Si+*) 

\ e,u iV(Si+) ^ ' 

where we have ignored coUisional de-excitation. The effects of collisional de-excitation are negligible when 
rie Her, where the critical density, ricr, is given by ricr = ^21/721- For T = 8000 K the critical density is 
ficr ~ 1600 cm~^. 

Using the integrated column density of Si II in component 1 yields (ne,i) = 0.2 cm~^. This result 
assumes T = 8000 K, and scales as T^/^. This value is a lower limit to the true electron density, ne,i. 
The Si+ absorption includes contributions from both H I and H II regions, while it is likely that the 
Si+* absorption comes primarily from only the densest gas of H Il-bearing regions. The column of Si+ 
associated with the ionized gas is overestimated, making Ue^i too small. The results of our photoionization 
modelling in ^.1| suggest that the column density of Si II arising in the H II region about /x Col is likely 
log7V(Si II) « 14.3. If we adopt this value for the column density of Si II associated with the ionized gas 
and assume that all of the Si II* arises in this H II, we can derive an upper limit for rie^i- Combining these 
two limits yields 

0.2 cm"^ <; 7ie_,; <^ 1.2 cm"^. 
The upper limit is larger than the values we've used for modelling the /x Col H II region, but the conclusions 



drawn in §4.1 are not sensitive to the ambient density. 



4-4-^- Electron Temperatures and Densities in the Neutral Gas 

We estimate the electron temperatures and densities of the primarily neutral material in the low- velocity 
components 1 and 2 by examining the ionization equilibrium of Ca and Mg in these clouds. The equation 
for ionization equilibrium of the ionic stages X'\ X*"*"^ of an element X may be written 

T{X')n{X') = a{X')nen{X'+^), (8) 

where a{X'^) is the recombination coefficient of to X^ and T{X^) the ionization rate of X\ The value 
a(X') is the sum of the radiative recombination coefficient, arad{X^), and the dielectronic recombination 
coefficient, adiiX'), and is a function of temperature. Using Equation (||) simultaneously for Ca^/Ca'^'^ 
and Mg°/Mg+ ionization equilibrium can yield estimates of both the electron temperatures and densities. 

For the recombination coefficients we adopt the fits suggested in the compilation 
of atomic data by D. Verner^. Namely for the Mg° recombination coefficients we use 
ad^(Mg°) = 4.49 X lO-'^T'^/^ exp(-5.01 x lO^/T) [1 + 2.1 x 10^^ gxp(-2.81 x lO^/T)] cm^ s'^ and 
arad{Mg°) = 1.4 X 10-13 (T/lO^y^'^^^ cm^ s'^ from ShuU & Van Steenberg (1982) and Aldrovandi & 
Pequignot (1973), respectively. We assume a photoionization rate r(Mg°) = 4.0 x lO^^^ s^^ (Frisch et 
al. 1990). For Ca we adopt the values a(Ca+) w a™d(Ca+) = 6.78 x lO'^^ (T/10^)~°'* cm^ s^^ from 
ShuU & Van Steenberg (1982)|. We assume r(Ca+) = 2.0 x lO-^^ g-i (pequignot & Aldrovandi 1986). 
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9 adi(Ca+) < a^ad(Ca+) for T < 10^ K. 
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Unfortunately, we do not measure the dominant ionization stage Ca+'^. To proceed we assume that the 
gas-phase abundances of Ca and Fe are approximately the same in the WNM towards /i Col (see Jenkins 
1987). The profiles of the Ca II and Fe II (e.g., A2249) transitions are very similar for the fi Col sightline, 
suggesting this is a reasonable assumption. Using this approximation we can substitute n(Fe~'') • {Ca/Fe}© 
for n(Ca++). 

Figure |l^ shows the relationship of and for Ca and Mg ionization equilibrium in component 
1. The solid lines represent the dependence of Ue on Te for each diagnostic given the best values of 
the iV(Ca+)/7V(Fe+) and 7V(Mg°)/iV(Mg+) ratios. The dotted lines give the la error limits based 
upon the sources of error given in § |3.2[ While there is a formal solution for very low temperatures, the 
high-temperature solution is more applicable to the gas being considered here. The results of this approach 
for components 1 and 2 are given in Table |l|. Also shown are the electron temperatures and densities for 
component 4 (see below). The treatment of Equation (||) neglects the charge-exchange reactions between 
neutral Mg and ionized H (Allan et al. 1988), but we have found including these reactions makes little 
difference to the final results. 

In principle one can derive the neutral hydrogen density in WNM clouds using observations of C I 
absorption from the ^Pq, ^Pi, and ^P2 levels. However, our marginal detections of the first excited ^Pi level 
and limits on the ^P2 column densities allow us only to place limits on the neutral hydrogen density. We find 
logiVio = log7V(C I 3Pi) -logAr(C I 3Po) < -0.37 and logiV2o = logiV(C I ^P2)-logN{C I ^Pi) < -0.56 
in component 1. Interpolating between the values given in the tabulation of Keenan (1989), and assuming 
T PS 6000 K and rie « 0.3 derived above, we estimate nn < 10 cm~'^ in component 1. 



5. IONIZATION AND DUST IN THE INTERMEDIATE- VELOCITY GAS 

5.1. Photoionized Gas at Intermediate Positive Velocities 

Components 3 and 4 along the sightline to /Lt Col are examples of low total column density intermediate- 
velocity clouds (IVCs). The hydrogen column density for component 4 is likely logAf(H) ^ 17.3 — 17.7 
(see below). With an absorption cross-section of 6.3 x 10^^^ cm^ for atomic hydrogen (Spitzer 1978), this 
implies an optical depth of order unity for this cloud at A912. Thus, the assumption that the gas-phase 
abundances of a cloud are accurately traced by the dominant stage of ionization in the neutral medium may 
break down for components 3 and 4. Further, the effects of ionization will differ from species to species, 
depending upon that ionic species' ionization potential and photoionization cross section, its recombination 
coefficient, including the effects dielectronic recombination, and the availability of charge- exchange reaction 
pathways. 

The species for which we have well-determined column densities for component 4 include N I, O I, 
Mg I and II, Si II and III, Fe II, and Al II. In Figure |l^ we plot the the values of [X/Si] for the species O, 
Mg, Si, Al, Mn, and Fe. The values [^/Si] for components 3 and 4, derived from the ratios of the usually 
dominant species, are also given in Table |^. The unusual nature of this component can easily be seen in 
Figure |l3|. Typically Si is moderately incorporated into grains in the diffuse ISM (e.g., in components 1 
and 2), while N has only moderately sub-solar gas phase abundances in the solar neighborhood (Meyer et 
al. (1997) find [N/H] « -0.1). Therefore in the WNM we expect to find [N/Si] > 0. In Figure |l| we see 
that [N/Si] derived from measurements of N I and Si II is sub-solar. At the same time we find that [Fe/Si] 
is almost indistinguishable from that seen in component 2 and similar to halo cloud values. We believe 
that the pattern of gas-phase column densities observed in component 4 is a result of photoionization of a 
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low column density cloud with somewhat typical halo- like gas-phase abundances. The uncertainties in this 
ionization, however, hinder our ability to make firm conclusions regarding the gas-phase abundances of the 
intermediate velocity clouds towards ^ Col. 

For a few species we can place quite stringent limits on the true gas-phase abundances. We observe the 
ionization states Si II, III, and IV and therefore have a very good measure of the total column density of Si 
in component 4, log A^(Si) = 12.85 ± 0.03, i.e., the contributions from Si III and IV are negligible. We can 
place similar constraints on the abundance of Al in this material. The 2a upper limit log A^(A1 III) <^ 11.2 
in component 4 implies —0.51 <^ [Al/Si] <j —0.25. If we assume Si is present in its solar system abundance, 
the column density of Si suggests the total hydrogen column density is logA^(H) « 17.35, whereas if its 
abundance is as low as [Si/H] w —0.3, i.e., having warm disk cloud-like abundances, then the total hydrogen 
column density of this cloud is logA^(H) w 17.65. 

Before discussing the effects of photoionization, we note that the observed gas-phase abundances 
are likely not due to coUisional ionization, nor is it likely these clouds have been ionized by the hard 
radiation from a strong shock. Recent works studying high-velocity material in the Galactic disk (Trapero 
et al. 1996) and gas associated with the Vela supernova remnant (Jenkins et al. 1998) have suggested 
these physical mechanisms are responsible for unusual abundances along the sightlines they studied. Two 
lines of evidence lead us to believe these are unlikely situations for the intermediate positive velocity gas 
towards Col. First, the observed 6-values for the IVCs towards Col imply very low temperatures. We 



derive a temperature for component 4 of T = 4, 000 ± 700 K in §5.3. Second, both of these works find 
increasing gas-phase abundances as a function of ionization potential (IP). In Figure |l^ we show the gas 
phase abundances of components 3 and 4 towards /i Col referenced to Si. The bottom plot shows these 
same data with the results from Jenkins et al. and Trapero et al. overlayed. We see very little trend with 
IP, except for a turnover at high IP with Si III. Trapero et al. find logiV(Si II)/A^(Si III) « —0.6, we find 
logiV(Si II)/Af(Si III) « -t-2.1 for component 4. 

Given the sub-solar abundance [Al/Si] derived above, it is clear that component 4 contains dust, which 
is affecting the pattern of column densties measured in this cloud. However, the abundances [0/Si] and 
[N/Si] derived from measurements of O I and N I reveal that photoionization is also playing a role. This 
conclusion is based on an examination of the relative ionization cross sections and recombination coefficients 
for the observed atomic and ionic species. The ionization equilibrium of a species in its simplest form is 
written in a manner similar to Equation (H). Ignoring for now charge exchange reactions, we write 

where r(X*) is a function of the ionization cross section and incident radiation field. Ions that have higher 
ratios of the ionization cross section cr^ to recombination coefficient a{T) are more readily ionized by a 
given photon of frequency i^. Sofia & Jenkins (1998) have discussed in detail the role of photoionization for 
Ar I, which has a large cross section to ionizing photons. They find that the observed deficit of Ar relative 
to H towards several stars is likely due to the preferential ionization of Ar I to II relative to the ionization 
of H. The results of SJ also show, when using the photoionization cross sections of Verner et al. (1996) and 
recombination coefficients of ShuU & Van Steenberg (1982), that O I and N I have significantly larger values 
r(A"*)/a(A'*) than most of the species studied here, while Al II, Mg II, and Si II were among the lowest. 
We suggest the deficits of N I and O I relative to Si II are the result of photoionization of this low column 
density cloud, similar to the deficit of Ar I relative to H I measured by SJ. The magnitude of the deficits is 
such that O would be ~ 50% ionized in component 4. We expect the ionization fractions of H to be very 
similar given the very strong charge exchange reactions between these elements (Field & Steigman 1971). 
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The fraction of C in the form 0+ is likely quite large given the large IP of this ion and the relative 
similarity of r(C^)/Q!(C^) to that of Si+ over a large range of energy. Sofia et al. (1997) have shown the 
gas-phase abundance of C is relatively constant at a value [C/H] « —0.4; if we assume this is also the case 
in component 4, then we expect [Si/H] w —0.3. However, the most important abundance measurement for 
this cloud is the value [Al/Si] derived above. Given our ability to account for the various stages of ionization 
of those elements, our measurements securely show this cloud contains dust grains. 

Component 3 may be even more prone to ionization effects. The column density of Si II in component 
3 is similar to that in component 4, though the column density of O I is lower by a significant amount, 
suggesting a higher degree of ionization {N{0~^~^)/N{0) « 0.7). The total column density of Si in this 
component is log A^(Si) = 12.83 ± 0.03, using the measured column densities of Si II and III and the upper 
limit for Si IV: log Ar(Si IV) < 11.3 (2a). The ratio log Af(Si II)/iV(Si III) = +0.95, while less than that for 
component 4, shows that the higher stages of ionization add relatively little to the abundance of Si. Again, 
we can accurately derive the gas-phase abundance of Al relative to Si given our upper limits to the column 
density of Al III. We find -0.31 < [Al/Si] < -0.11 for component 3. Thus the values [Al/Si] may be slightly 
different in components 3 and 4. Given our 3c7 upper limits to the column densities of S II and III in this 
component, [Si/S]^ —0.15, or —0.3 if we assume the amount of S III is negligible. This suggests component 
3 could be similar in its gas-phase composition to component 2. 

5.2. CoUisionally Ionized Gas at Intermediate Negative Velocities 

Component 5 likely samples a low-column density region of warm, coUisionally ionized gas. For the 
species that have well-determined 6- values, i.e., C II, Mg II, and Si III, we derive upper limits to the 
temperature of T < 36, 000 K, < 82, 000 K, and < 95, 000 K, respectively. Unfortunately, our 6-values are 

not well enough constrained to disentangle the thermal from non-thermal broadening. For C II and Mg II, 
ions with a factor of two difference in atomic mass, to have such similar 6-values, non-thermal broadening 
must play a significant, if not dominant, role in this component. 

For component 5 we find logA''(Si II)/7V(Si III) = —0.67, similar to the value observed towards 23 Ori 
by Trapero et al. (1996). This value is also similar to a gas in coUisional ionization equilibrium having a 
temperature slightly less than T w 25, 000 K (Sutherland & Dopita 1993), which consistent with our upper 
limit of 36,000 K. If Si is not depleted in this gas and all of the Si is found in the Si"*" and Si+^, then the 
column density of hydrogen (H° -|- H+) in this component should be log A'^(H) ^ 16.8. We quote a lower 
limit since some depletion is likely. 

Photoionization models using stellar input spectra or models for the interstellar radiation field are 
unable to match the observed ratio of Si II to Si III (e.g., Howk & Savage 1999) imless the gas is directly 
exposed to a hard radiation field. It is thus unlikely that this component represents H II region gas or gas 
associated with the WIM, even if it is made up of several overlapping clouds. Howk & Savage (1999), for 
example, find log7V(Si II)/A/'(Si III) ^ 0.0 for H II regions surrounding stars having T^ff ^ 39, 000 K. 

If component 5 is indeed in coUisional ionization equilibrium with T Ri 25, 000 K, the fraction of 

C in stages other than C"*" should be minimal given the high IP of C"*" (Sutherland & Dopita 1993). 
Using log[iV(Si II) -|- iV(Si III)] = 12.37 ± 0.03 as the column density of Si in this component, we find 
[C/Si] = 0.0 ±0.1. Various authors have found gas-phase abundances of [C/H] w —0.4 in the local ISM 
(referenced to solar system values) with very little intrinsic spread (Sofia et al. 1997; Cardelli et al. 1993). 
Thus, if we assume this value also holds in component 5, our derived [C/Si] suggests [Si/H] similar to that 
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found in component 1. 

We also note the presence of an absorbing component at wlsr = —49.5 ± 0.3 km s^^ visible 
only in Si III. This component may be similar to component 5. Component fitting analysis yields: 
logA^(Si m) = 11.62 ± 0.03, and 6 = 5.9 ± 1.0 km s"\ constraining the temperature T ^ 59,000. 

The relative ionization levels of component 5 are similar to the clouds observed towards 23 Ori by 
Trapero et al. (1996), as evidenced by the relative column densities of Si II and III. There are two main 
observational differences between the collisionally-ionized material along the /i Col and 23 Ori sightlines, 
which may or may not imply true physical differences between these clouds. First, the collisionally-ionized 
gas towards 23 Ori is at high velocities relative to the LSR of ?;lsr ~ —100 to —125 km s~^. However, while 
component 5 does not have as extreme a velocity as the ionized material seen towards 23 Ori, its velocity 
is inconsistent with models of Galactic rotation in this direction. This alone, however, is not extremely 
unusual given that components 2, 3, and 4 towards /i Col are also at odds with simple models of Galactic 
rotation. These clouds make up only ~ 10% of the total column density along the fi Col sightline. 

A second difference between the material making up our component 5 and the HVCs observed by 
Trapero et al. is the observed component structure of the gas. Trapero et al. report the existence of four 
well-separated components in the HVC seen towards 23 Ori, each of which is characterized by a 6-value 
between 2 and 3 km s~^. Thus they derive temperature limits T <^ 12, 000 K for each of their components. 
These authors suggest the gas is not in collisional ionization equilibrium, but rather that the ionization 
state of the gas is "frozen in" with relative ionization levels more representative of T 25,000 K. The gas 
temperature suggested by comparing the amounts of Si II and III with ionization equilibrium calculations 
for component 5 ('^ 25,000 K) is consistent with the 6- values derived from our component-fitting analysis 
(suggesting T < 36,000 K). We must caution, however, that although our observations of component 5 
are made with the echelle-mode gratings of the GHRS, with a resolution of « 3.5 km s~^, there may 
be unresolved component structure present in the profile. Though this component is reasonably well 
fit by a model cloud represented by a single gaussian, this does not exclude the presence of multiple 
components within this velocity range. Indeed, the similarity of the C II and Mg II 6-values suggests that 
the temperature of the gas in this component is much less than the upper limit derived from our the C II 
fe- value. 

Thus the intermediate negative velocity gas along the fi Col sightline appears to be an example of low 
column density, coUisionally ionized material. The two clouds at wlsr ~ — 29 and —49 km may be 
similar to the HVCs observed towards the disk stars 23 Ori and r CMa. Though the negative velocity IVCs 
towards /i Col have slightly lower column densities and appear at lower velocities, the ionization state of 
these IVCs appear to be similar to the HVCs presented by Trapero et al. (1996). 

5.3. Physical Conditions in the Intermediate- Velocity Gas 

Our component fitting results for component 4 include good 6-values for species over a wide range of 
atomic mass, from N I to Fe II. Figure |l^ shows a plot of 6^ vs. I /A, where A is the ion mass in atomic 
mass units, for component 4. The best linear fit to these data is also overplotted. Theoretically this fit is a 
sum of the expected broadening from thermal and non-thermal sources such that: 



(10) 
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where toh is the mass of hydrogen and Vnt is the non-thermal velocity dispersion. We determine a 
temperature T = 4, 000 ± 700 K with a non-thermal velocity component Vnt = 0.6 ± 0.4 km s~^ for 
component 4. While the fit to the data is quite good, this temperature is slightly lower than typically 
found in diffuse interstellar clouds (e.g., T « 6,000 K by Spitzer & Fitzpatrick 1993). If we use only the 
measurements for Fe II, Si II, and Mg II in this fit, we derive T = 4, 300 ± 700, in agreement with the values 
derived when including the less certain O I and N I &- values. 

Spitzer & Fitzpatrick (1993) have discussed the population of the upper fine-structure level of the C+ 
ground state via electron collisions. Their Equation (7) relates the value of rif, in diffuse clouds to the 



column of C+* in the same way we have approached the excitation of the Si^ lines in §4.4. The extreme 
strength of the C+ line at A1334 makes the column density of the lower-lying level uncertain. Adopting the 
atomic constants used by Spitzer & Fitzpatrick (1993), namely: A21 = 2.29 x 10~^ s~^ (Mendoza 1983) and 
= 2.90 (Osterbrock 1989), we find 

(11) 

We have assumed that the column densities may be used in place of the particle densities. Using 
Equation ( pT| ) with the temperature T = 4, 000 ± 700 K derived from our component fitting, and assuming 
log A^(C+) = 13.75 ± 0.1, gives = 0.47 ± 0.14 cm~^. Analysis of the Mg°/Mg+ ionization cquihbrium 
using Equation and this temperature yields rig = 0.64 ± 0.14 cm~^, in surprising agreement with the 
density derived from C'*'* . 

The cooling per nucleon via the C II ^^3/2 — P1/2 radiative transition at A158 /xm can be written 

l, = hi^i2Ai2N{CU*)/N(H), (12) 

where hi'12 = 1.26 x 10^^ ergs is the energy of the emitted photon, and A21 = 2.29 x 10^^ s^^ is the 
Einstein A- value for the transition (Mendoza 1983). If log A^(H) « 17.3 — 17.7 in this component, our C II* 
column density yields /c = (3 — 4) x 10~^^ ergs s^^ nucleon^^. This value is similar to that derived for 
sightlines in the Galactic disk (Pottasch et al. 1979; Gry et al. 1992), and higher than observed along most 
halo sightlines (Savage & Sembach 1996b). 

The charge exchange reaction between 0+ and H° is large enough to effectively keep the ionization 
fractions of O and H coupled. If we assume that O and Si in components 3 and 4 have the same intrinsic 
gas-phase abundances (Si in low-density regions may be slightly more abundant), we can use the observed 
deficit of O I to estimate the ionization fraction of H in these clouds. Using the values [0/Si] in Table ^ 
and again using Xe = n^/n-n ~ A^(H+)/A^(H), we find 0:3 w 0.7 and k, 0.5. If we assume ~ 0.5 cm~^ 
in component 4, we estimate nu « 1.0 cm~'^, assuming He is mostly neutral. The corresponding thermal 
pressure is P/fc « 6,000 K cm~^, similar to other pressure estimates for diffuse gas. The assumption of an 
intrinsic value [0/Si] = 0.0 likely makes the estimated ionization fractions upper limits and the densities 
and pressures lower limits. 

If we assume x ~ 0.5 and rif. w 0.5 cm^'^ for component 4 and adopt the column density for H derived 
above (logiV(H) « 17.3 — 17.7), the thickness of the absorbing region is ^ 0.04 — 0.1 pc. For the HVC 
towards r CMa, Trapero et al. (1996) suggest the absorbing region is likely ^ 0.065 pc. 



6. DISCUSSION 
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6.1. Interstellar Reference Abundances 

An important result of this paper is our conclusion that the relative abundances of S, P, and 
Zn are consistent with solar system relative abundances, i.e., [Zn/S] « [P/S] ~ 0.0. Furthermore, 
[S/H] +0.08 ± 0.02, [P/H] = +0.05 ± 0.02, and [Zn/H] = +0.06 ± 0.08 using the integrated sighthne 
column densities with no correction for the presence of ionized gas along the sightline. There is no evidence 
for saturation in our S II observations of the weaker two lines (AA1250.534 and 1253.811; see Figure 



and the strongest line shows only moderate saturation effects. In §4.1 we showed that the contribution of 
ionized gas in an H II region to these species was in the range 0.04 to 0.05 dex. An additional correction 
of ~ —0.04 dex may need to be applied, in particular, to the S II measurements to correct for the presence 
of partially-ionized material in the predominantly neutral H I-bearing clouds. Thus relative to H, the 
abundances of S, P, and Zn are consistent with solar (within the errors). 

The conclusion that S, P, and Zn relative to H are consistent with their solar system abundances 
(or greater) is significant because we believe we understand the effects of ionized gas on these and other 
species along this sightline. Abundance measurements of [S/H] derived from S II and H I are particularly 
susceptible to contamination from H+-containing regions. In ionized regions, which are not included in the 
derived H I reference column density, the dominant ion of S can be S+ given its high ionization potential 
(23.3 eV; see Howk & Savage 1999). Thus the derived gas-phase abundance [S/H] may be unreliable if 
significant amounts of H+ are present along a sightline. Comparing our quality measurements of H I with 
Zn II and P II, which have lower ionization potentials (18.0 and 19.7 eV, respectively) than S II, allows us 
to firmly conclude that we are seeing gas with solar system abundances in S, Zn, and P. 

The derived solar system relative abundances of the gas towards ^ Col are not tracing initially 
sub-solar metallicity material contaminated by a SN or other enrichment event. The measurement of S, 
an a-process element, Zn, a tracer of Fe-peak elements, and P, an odd-Z element produced primarily in 
the hydrostatic O- and Ne-burning shells of massive stars (Timmes et al. 1995), shows no evidence for 
a different nucleosynthetic signature in the gas towards /i Col than that responsible for producing the 
observed metals in the solar system. 

Recent studies by Sembach et al. (1995) and Roth & Blades (1995) have shown the abundance of Zn in 
the ISM is a function of the sightline value of log A(H I) as well as the fractional abundance of molecular 
hydrogen, /(H2) = 2A^(H2)/[2A^(H2) + A(H I)]. This dependence on sightline parameters such as column 
density, average sightline neutral density, and molecular content may be interpreted as evidence for the 
incorporation of Zn into grains in the ISM (e.g., Jenkins 1987). The n Col sightline samples low density 
WNM gas {{iim) = -/V(H l)/d ^ 0.06 cm~'^). Therefore it would not be surprising if this sightline were to 
show higher gas-phase abundances of Zn than a sample including high density sightlines. However, both 
Sembach et al. (1995) and Roth & Blades (1995) find values ([Zn/H]) w -0.2 in their sample of Milky Way 
disk and halo stars with little molecular gas [log/(H2) ^ —3]. This is significantly lower than the observed 
[Zn/H] towards fi Col, which has log/(H2) ~ —4.35 according to the results of SY. The H I column density 
towards /i Col is significantly lower than any in the Sembach et al. sample and lower than all but one of 
the stars studied by Roth & Blades. The few stars from the latter sample with average sightline neutral 
densities similar to the sightline to /x Col still show [Zn/H] ~ —0.1, though the errors are likely of order 0.08 
dex. 

Quality, high- resolution measurements of [S/H] are somewhat more sparse in the literature, in part 
because the transitions of S II are stronger than those of Zn II and therefore saturate along even relatively 
low column density sightlines. We have collected from the literature the following quality integrated sightline 
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measurements of [S/H] using the S II transitions: [S/H] = 0.00 ± 0.09 towards 7^ Vel (Fitzpatrick & Spitzer 
1994); [S/H] = -0.05 ± 0.04 towards HD 93521 (Spitzer & Fitzpatrick 1993); and [S/H] = -0.09 ± 0.09 
towards HD 215733 (Fitzpatrick & Spitzer 1997). These values were derived assuming H I column densities 
from Diplas & Savage (1994), except for the 7^ Vel [S/H] value, which uses an H I column density that is the 
weighted average of the Diplas & Savage (1994) measurement with those of Bohlin et al. (1978) and York & 
Rogerson (1976). All of these measurements are roughly consistent with a solar system abundance of sulfur. 
The number of measurements is too few and the errors on many of the measurements are too large to 
provide good estimates of the intrinsic scatter in the S abundance. We do note that the high-quality fj, Col 
and HD 93521 measurements are inconsistent with one another at 2(t significance. However, we have not 
accounted for the effects of ionized gas. We have seen that towards /j. Col the contribution from ionized gas 
to the S II measurements can be of order 0.05-0.08 dex. The values given above should be viewed as upper 
limits to the true gas-phase abundance of sulfur. Disentangling the influence of ionized gas on the sulfur 
abundance for many sightlines, where the fraction of ionized gas may vary considerably, may be non-trivial. 

Our measurements of solar-like abundance ratios along the /i Col sightline argue for a solar reference 
abundance along at least one sightline. Comparing the [Zn/H] abundance towards /i Col with [Zn/H] along 
the low-density sightlines studied by Roth & Blades (1995) and Sembach et al. (1995) suggests possible 
abundance variations in the low-density ISM. Our measurements of solar system relative abundances along 
this sightline are in striking contrast to the results of Meyer et al. (1997, 1998) and Cardelli & Meyer (1997) 
regarding the constancy of (sub-solar) 0/H, N/H, and Kr/H abundances in the ISM within <^ 1.5 kpc of 
the sun. Our data are inconsistent with [Zn/H] = —0.2 or [P/H] = —0.2 (for comparison with the results of 
Sembach et al. and Roth & Blades, as well as the Meyer et al. [0/H] results) at the 2.5a and lOcr levels, 
respectively, after subtracting off the expected contribution from ionized gas (0.04 to 0.05 dex). 

There has been growing evidence in the literature to support the choice of a sub-solar reference 
abundance in the local ISM (Meyer et al. 1998; Mathis 1996; Snow & Witt 1995, 1996). The arguments 
for sub-solar abundances appear to be supported by observations of local B star abundances (e.g., Gies & 
Lambert 1992; Kilian-Montenbruck et al. 1994) and H II region abundances (e.g., AfFlerbach, ChurchweU, 
& Werner 1996; Simpson et al. 1995; Osterbrock, Tran, & Veilleux 1992; Shaver et al. 1983). However 
our results suggest the issue is not resolved. There may be intrinsic abundance variations in the ISM 
near the sun. For many sightlines the solar system abundances may actually be the appropriate reference 
abundance system. Some caution is warranted given that we have identified only one sightline with solar 
system abundances. However, the fi Col sightline is exceptionally well observed, and our understanding of 
the ionization state of the gas along this path through the ISM is far better than for most others. The 
interpretation of the measured 0/H and N/H ratios are simplified by the strong charge exchange reactions 
between 0° and H°, and to a lesser extent N° and H°. This keeps the relative ionization fractions of these 
species locked to that of H. Though our measurements are of species less strongly tied to hydrogen, the fact 
that P and Zn appear in their solar system ratios relative to H, and that our photoionization modelling can 
account for the observed super-solar abundance of S/H, gives us confidence that we are indeed seeing solar 
system abundances along this sightline. 

The identification of the proper reference abimdance for the local ISM, or for an individual sightline, 
is likely not an easy task. Furthermore, the observed abundances for the B-star and H Il-region reference 
systems are less well constrained at this time than for that of the sun. Indeed, the spread in the B-star 
abundances (0.2-0.7 dex) seen by Kilian-Montenbruck et al. (1994) may be evidence that there are intrinsic 
variations in the local reference abundance. However, these authors see 0.2 to 0.7 dex abundance spreads 
for stars within the same clusters, possibly suggesting very small scale abundance inhomogeneities in the 
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material from which the stars were formed or inhomogeneities introduced by the star formation process 
itself. At this point we caution that it is not clear that either solar system abundances or some constant 
fraction thereof should be adopted out of hand as the cosmic reference abundance. 



6.2. Implications for High-Redshift Absorption Systems 

Our work has important implications for the study of high-redshift quasar absorption line systems, 
particularly the damped Lya systems (DLAs). Among the most important aspects of this paper for 
interpreting absorption line spectroscopy of higher-redshift gas are the solar system abundances of 



undepleted elements as discussed in §5.1, and in particular the measurement [P/Zn] « 0.0, and the 



application of new Ni II oscillator strengths to the study of the gaseous medium of the Galaxy. 

The new high-quality laboratory measurements of Ni II absolute oscillator strengths imply that 
previously derived Ni column densities should be increased by 0.3 — 0.4 dex, depending on the source 
of the /-values. In particular, column densities derived using the /-values of Zsargo & Federman (1998) 
should be increased by -1-0.272 dex, although there may still be uncertainties in the oscillator strengths of 



the lower-wavelength transitions (e.g., AA1317 and 1370; see §3.2). Nickel can be an extremely important 
element for studies of DLAs because the Ni/Fe ratio is set by nuclear statistical equilibrium in the sites 
of Ni and Fe production. Thus the ratio of Ni/Fe (dust-|-gas) is expected to be the same no matter 
the nucleosynthetic history of the material. Observations of disk and halo stars over a large range of 
metallicities confirm this expectation (Edvardsson et al. 1993; Gratton & Sneden 1988, 1991; see also 
discussion in Lu et al. 1996), although there has been disputed evidence for a slight increase in [Ni/Fe] 
at very low metallicities (see discussion in Wheeler, Sneden, & Truran 1989 and McWilliam 1997 and 
references therein). Edvardsson et al. (1993) derive [(Ni/Fe)] = +0.02 ±0.05 for their sample of 189 galactic 
disk G and F dwarfs, which cover a range of abundances —1.1 [Fe/H] -fO.25. Gratton & Sneden (1991) 
similarly find [(Ni/Fe)] = -0.05 ± 0.08 for their sample of 22 metal-poor stars ([Fe/H] ^ -0.8). 

Measurements of the gas-phase [Ni/H] in disk and halo clouds by numerous authors have suggested 
[Ni/Fe] < 0.0 (see Savage & Sembach 1996a and references therein) when using the Ni II /-values from 
Morton (1991). Several groups have found sub-solar Ni/Fe ratios in DLAs as well (e.g., Lu et al. 1996; 
Prochaska & Wolfe 1999; Pettini et al. 1999). This depletion of Ni relative to Fe in DLAs has been 
interpreted to imply the existence of dust grains in these systems, given that one expects [Ni/Fe] « in the 
absence of dust. Lu et al. (1996), however, pointed out the possible uncertainties in the /-values of the Ni II 
transitions and adopted to reject this line of evidence for presence of dust given the uncertainties in the 
atomic data. The recent Fedchak & Lawler measurements of the Ni II resonance oscillator strengths revise 
all of the previous DLAs Ni abundance measurements upwards by ^ 0.3 dex. This is typically enough to 
bring the Ni/Fe ratio near the solar system abundance in DLAs. Combining the 20 DLAs with detected 
Fe and Ni absorption from the observations of Lu et al. (1996) and Prochaska & Wolfe (1999) we find 
[(Ni/Fe)] = 0.15 ± 0.25 in these DLAs (typical measurement errors are ~ 0.1 dex). The absorbing systems 
in this combined dataset cover a range in redshift 1.8 ^ z ^ 3.9; the [Ni/Fe] measurements seem not to be 
a function of redshift. 

The revisions in the Ni II oscillator strengths imply that [Ni/Fe] measurements no longer offer a 
straightforward indicator of the presence of dust. We calculate a revised value [Ni/Fc] = +0.04 ± 0.05 
for cloud A along the C Oph sightline (using both revised Ni II and Fe II oscillator strengths). Coupled 
with our own measurements of [Ni/Fe] — +0.02 ± 0.06 towards ^ Col (see Table |[), this implies the ratio 
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Ni/Fe can be consistent with solar system abundances, even in the presence of dust. Many of the existing 
claims for dust in DLAs that use sub-solar Ni/Fe ratios as evidence should be reconsidered given the new 
absolute /-value scale for the AA1710 and 1741 transitions, and the suggested relative scale for many other 
transitions. This result once again emphasizes the need for high-quality oscillator strength measurements. 

Another possibly important aspect of our work for studies of high-z gas is the observation [P/Zn] « 
along the fi Col sightline. The elements S and Zn have often been used to search for enhancements of 
a-process elements over Fe-peak elements in the distant universe, particularly because they are not expected 
to be incorporated into dust grains in large amounts. These comparisons offer clues to the nucleosynthetic 
history of the gas, which can be compared with that of low-metallicity Galactic stars. Unfortunately, the 
strength of the S II transitions relative to the Zn II lines often means that the S lines are saturated where 
the Zn lines are strong enough to be detectable. The P II A1152 transition is expected to have a strength 
closer to the Zn II (e.g., see Table ||) and might therefore be better choice for comparing abundances of an 
elemental produced by type II SNe (P) with one tracing Fe-peak elements (Zn). 

The interpretation of the P/Zn ratio is not as straightforward as that of S/Zn, however. First, 
a-elements, such as S, are known to show enhancements relative to Fe-peak elements, traced by Zn, in 
low-metallicity Galactic halo stars (e.g., Wheeler et al. 1989 and references therein). This is understood 
as a result of the time delay between the explosions of massive stars as type II SNe (which produce the 
a-process elements) and the explosions of white dwarfs in binary systems as type la SNe (which in turn 
produce much of the Fe-peak elements). Unfortunately the emperical information about the chemical 
evolution of P is not as well constrained given the lack of abundance measurements of P in Galactic stars 
(see Timmes, Woosley, & Weaver 1995). The production of '^^P primarily occurs in the hydrostatic O- and 
Ne-burning shells of pre-supernova massive stars (Woosley & Weaver 1995; Trimble 1991). Unlike many 
elements, the amount of '^"'^P produced in a type II SN is relatively insensitive to the exact mechanism of 
the explosion (Timmes et al. 1995). Timmes et al., in their detailed chemical evolution calculations, find 
the abundance [P/Fe] should be significantly sub-solar in low metallicity environments ([P/Fe] ^ —0.8), and 
rise to super-solar values ([P/Fc] +0.15) as the metallicity [Fe/H] increases. Therefore, the yield of P in 
type II SN explosions is expected to be a function of the initial metallicity. This is similar to the expected 
behavior of [N/Fe] if the production rate of nitrogen depends on the initial CNO abundances, i.e., if the 
synthesis of N has a secondary origin. The measurement of [P/Zn] in high-rcdshift DLAs thus probes the 
chemical evolution of an element with a secondary origin in massive stars and an Fe-peak element produced 
primarily through low-mass stars. 

The validity of conclusions drawn from P and Zn comparisons in DLAs depends on the assumption that 
P and Zn are undepleted in the gas, or have similar levels of depletion. Our measurements suggest that the 
former is true along a low-density WNM sightline in the Galactic disk. The abundance of P in the diffuse 
ISM of the galaxy has not been well characterized. Dufton, Keenan, & Hibbert (1986) used Copernicus data 
to study the abundance of P; they concluded that the abundance of P in low-density clouds was consistent 
with solar. However, there were only a handful of P II measurements in their sample of warm, low-density 
sightlines. Savage et al. (1992) find [P/Zn] « —0.2 in the warm cloud A towards Q Oph (using our adopted 
oscillator strengths), which otherwise shows abundances similar to component 1 along the /i Col sightline. 
While the /i Col data suggest that P and Zn trace each other well in the Galactic ISM, the contradictory 
results for a similar cloud towards C, Oph suggest a better sample of P measurements is needed to constrain 
the relative behavior of P and Zn in the Milky Way. Measurements of the P to Zn abundance in the ISM 
of the Magellanic Clouds and other local low-metallicity environments may also be helpful in this regard. 
Without this fundamental knowledge of the properties of low-redshift gas, the usefulness of P and Zn as 
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tracers of high-redshift chemical evolution may be limited. 



6.3. Composition of Dust in the Diffuse ISM 



In Table 11 we have compiled data for the dust-phase abundances of Mg, Si, and Fe in well-studied 
diffuse interstellar clouds. The data are taken from the works of Fitzpatrick & Spitzer (1997), Savage et al. 
(1992), Sofia et al. (1994), and Spitzer & Fitzpatrick (1993, 1995) as well as this work. References to the 
abundance measurements are given in the table. We only include clouds for which quality measurements 
of Mg, Si, and Fe exist. The data have been adjusted to reflect our choice of Mg II oscillator strengths 
(Fitzpatrick 1997), and reference abundance (see Savage & Sembach 1996a). The table has two sections: 
one showing the dust-phase abundances of several "warm disk" clouds for comparison with the dust-phase 
abundances in component 1, and the other containing information on "halo" clouds for comparison with 
component 2 (see Sembach & Savage 1996). Also given are the weighted mean of each of these samples and 
the standard error of the mean. The sample presented here does not represent an unbiased selection, but 
the averages are useful for judging the significance of variations in the dust-phase abundances within these 
clouds. 

Our selection of clouds for each group (halo or disk) is based upon the Sembach & Savage (1996) 
categorizations and the similarity of sightline environments. As such we have chosen clouds with similar 
properties for each category. The small spread in (Mg/H)^, (Si/H)^, and (Fe/H)^^ in the warm disk clouds 
listed in Table |l^ is quite interesting. The dust phase abundances of Mg, Si, and Fe show a spread of 
< ±10% from the mean values, with the exception of the (Fe/H)^ measurement for HD 18100. This 
sightline, for which we have used the sightline-integrated values of [X/Zn] in deriving the dust-phase 
abundances (Savage & Sembach 1996b), likely contains a larger fraction of low-density material than the 
others categorized as warm disk clouds. Even given the large difference in (Fe/H)^ for HD 18100 [and 
correspondingly large excursion of (Fe/Si)d from the mean], the value of [(Mg-|-Fe)/Si]d is still only ^ 10% 
from the mean value. Thus it would appear that a population of warm clouds in the Galactic disk can be 
identified with very similar dust composition. Component 1 towards /i Col seems to lie on the low end of 
the (Mg/H)d and (Si/H)c( values, which is consistent with the slightly higher [Mg/S] and [Si/S] values in 
this blend compared with cloud A towards C Oph. The dust-phase abundance of Fe is close to that observed 
in other warm disk clouds. The dust-phase abundances associated with the disk clouds in Table |ll| may 
be representative of grains in the WNM of the disk, particularly given the relatively small variation in the 
abundances. 

The data given in Table ^ suggest there is a quite small spread in the dust-phase abundances for 
clouds having "halo" -like abundances, although the errors in the available data are larger than for typical 
disk clouds. The dust-phase abundances of Fe, Mg, and Si relative to H are consistently lower in the halo 
clouds than the disk clouds of Table ^ Again, /i Col seems to lie near the lower end of the dust-phase 
abundances. This is consistent with the findings of SSC. In Figure 6 of Sembach & Savage (1996), which 
shows gas-phase abundances for many sightlines, the SSC results for component 2 along the fj. Col sightline 
typically show the highest abundance [A"/!!] for a given element. Though the errors are large, it does not 
appear that the ratios of Si, Mg, and Fe inclusions in the grains associated with component 2 are different 
than for the other halo clouds listed in Table 0. This cloud is therefore characterized by a greater return 
of elements to the gas-phase than the other halo clouds with no fundamental change in the dust-phase 
composition. 
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Table shows that the dust-phase abundances of both typical WNM disk clouds and halo clouds are 
inconsistent with the sole constituent of Galactic dust being pure olivine or pyroxene grains. Fitzpatrick 
(1997) has also pointed out that the gas-phase abundances from works of Fitzpatrick & Spitzer (1997) 
and Spitzer & Fitzpatrick (1993, 1995) are consistent with non-silicate inclusions of Fe and/or Mg into 
grains, possibly in the form of oxides, if adopting the Mg II /-values he derives in that work. It has been 
suggested that the dust-phase abundances in the halo clouds are representative of the composition of the 
resilient cores of grains, and in the warm disk clouds of the cores plus mantles (see Savage & Sembach 
1996a). If this is the case then the data in Table p] suggest the mantles contain Si, Fe, and Mg in the ratio 
Si:Fe:Mg « 7:5:4. Therefore, [(Mg-|-Fe) /Si] mant/e ~ 1-4, which is consistent with a mixture of the common 
silicates mentioned above. This is slightly larger than the value [(Mg-|-Fe)/Si]ma„t;e ~ 0.94 derived by 
Savage & Sembach (1996a). This difference is primarily due to the assumed Mg II oscillator strengths. The 
data in Table ^ imply the postulated resilient grain cores must include a greater fraction of non-silicate Fe 
and Mg inclusions than the more fragile mantles. 

The fact that the different oscillator strengths used for the Mg II doublet near A1240 variously imply 
no Mg depletion (e.g., Spitzer & Fitzpatrick 1993, 1995), moderate Mg depletion similar to that of Si (this 
work), and relatively large Mg depletion similar to that of Mn (e.g., Sembach & Savage 1996) shows the 
need for an an accurate laboratory measurements of the /-values for these important transitions. 



7. SUMMARY 

We present absorption line profiles from the GHRS for many atomic and ionic species along the line of 
sight to the nearby star fj, Col. The majority of these data were taken with the echelle-mode resolution of 
w 3.5 km s^^, which we analyze with a combination of the apparent column density method and component 
fitting techniques. This dataset is the most extensive interstellar absorption line database for studying the 
WNM of the Galaxy. The principal results of this study are: 

1. The presence of ionized gas near wlsr ~ km s^^ is suggested by the presence of S III, Al III, Fe III, 
Si III and IV, and Si II* . Copernicus observations of N II and N II* and WHAM spectra of Ha emisson 
also imply ionized gas at these velocities. Based upon CLOUDY photoionization models and our 
measurements of the doubly ionized species, we estimate the contribution from gas associated with 
ionized hydrogen to the column densities of singly ionized species used as abundance tracers in the 
neutral gas. We find that the relative contribution from ionized gas to these tracers is minimal. 

2. We find no evidence for sub-solar gas-phase abundances of Zn or P relative to S over the velocity range 
for which we are able to measure these species, i.e., [Zn/S]w[P/S]sa 0.0. The gas-phase abundances 
[Zn/H] and [P/H] are also solar within the errors. The [S/H] abundance derived from S II and H I 
is super-solar by -fO.08 dex, though our photoionization model, which corrects for S^ in ionized gas 
reduces this by —0.05 to —0.08 dex. This suggests the proper reference abundances for investigating 
the gas along the sightline to /i Col are the solar system abundances (e.g., Anders & Grevesse 1989). 

3. Adopting the /-values derived for Mg II AA1239.925 and 1240.395 by Fitzpatrick (1997), we find 
the abundance patterns of Mg and Si in the WNM towards /i Col as a function of velocity are very 
similar. This result is in accord with the trends noted by Fitzpatrick (1997). 

4. The gas-phase abundances for the warm neutral absorbing complex centered at wlsr = +3 km s^^ 
are similar to those found in the warm cloud towards C Oph and other warm disk clouds (Savage & 



- 33 - 



Sembach 1996a). Table | gives our derived gas-phase abundances which are shown in Figures and 
|lT| . This level of gas-phase abundances seems to be typical of WNM clouds in the diffuse ISM of the 
Galactic plane. 

5. Component 2 at wlsr = -1-20.1 km s~^ shows higher gas-phase abundances of the refractory elements 
compared to component 1 , suggesting grain processing has played a part in the evolution of this cloud. 
Comparison of the gas-phase abundances of components 2 to 1 suggests the liberation from the the 
solid phase (55 ± 20)% of Si, (50 ± 4)% of Fe, (23 ± 8)% of Ti, (60 ± 9)% of Cr, (26 ± 5)% of Mn, and 
(32 ± 4)% of Ni. We rule out enrichment from Type la SNe as the cause of these enhanced gas-phase 
abundances. 

6. The values of the dust-phase abundances of Mg and Fe relative to Si derived for the low- velocity 
material towards /i Col, [(Mg-|-Fe)/Si]d ~ 2.7 — 3.3, assuming solar relative abundances, are 
inconsistent with incorporation of Fe and/or Mg solely into olivine and pyroxene-type silicates, which 
predict [(Mg-|-Fe)/Si]£; = 1 — 2. Oxides are a likely component of the Mg- and Fe-bearing grains in 
the neutral ISM towards this star. This result is even more pronounced if one adopts B-star reference 
abundances. 

7. The low-velocity gas (components 1 and 2) along this sightline is characterized by Tg w 6, 000 — 7, 000 
K and rtg « 0.3 cm^'^, derived from ionization equilibrium of Mg and Ca. The temperature is typical 
for warm diffuse interstellar clouds of the Galactic disk and low halo (e.g., Spitzer & Fitzpatrick 
1993). In the densest ionized gas regions at low velocities we find 0.2 < rtg ^ 1.2 cm~'^ from excitation 
equilibrium analysis of the ^Pzii and fine structure levels of Si+. This ionized gas may trace an 
H II region about /i Col. 

8. The relative column densities of atomic and ionic species in the IVCs towards /i Col, components 3 and 
4 at v\sK = -1-31.0 and -1-41.2 km s~^, respectively, are dominated by the effects of (photo)ionization 
and elemental incorporation into grains. Our conclusions regarding these clouds are limited by the 
uncertain ionization state of the gas. However, we can bracket —0.31 <^ [Al/Si] <^ —0.11 in component 
3 and —0.51 <^ [Al/Si] <; —0.26 in component 4. These measurements are based upon several stages 
of ionization and should therefore be quite robust. Component 5 at wlsr = ^30 km s~^ traces 
coUisionally ionized gas and is only measured in five species. 

9. From our component fitting results for component 4 we derive T = 4, 000 ± 700 K. Using coUisional 
excitation equilibrium and the column densities of C II in the and fine structure levels, we 
derive rie = 0.47 ± 0.14 cm^"^. From ionization equilibrium of Mg I and II we derive rie — 0.64 ± 0.14 
cm~'^. The cooling in the C II A158 /zm line via radiative decay from the ^P3/2 level is likely in the 
range (3 - 4) x 10"^^ ergs s~^ nucleon"^. Component 5 has T <^ 36, 000 K. The ratio of Si II to III is 
similar to that expected for gas in coUisional ionization equilibrium at T « 25, 000 K. 
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A. DERIVATION OF iV(Hl) TOWARDS COL 

In this appendix we derive the column density of neutral hydrogen along the sightline towards /i Col. 
The H I column density for this sightline has been determined previously by Diplas & Savage (1994) 
and Bohlin et al. (1978) from lUE and Copernicus satellite spectra. However, the high signal to noise 
observations taken of the region near Lyman-a (Lya) with the intermediate resolution gratings of the 
GHRS, coupled with their exceptional scattered light properties, make a new determination of Af(H I) a 
useful endeavor given our final goal of deriving very accurate gas-phase abundances along this sightline. 

To derive the H I column density towards /i Col we use the continuum reconstruction method (Bohlin 
1975; Diplas & Savage 1994; Sofia et al. 1994). The optical depth of an interstellar absorption line as a 
function of wavelength, t\, is (Spitzer 1978): 

2 

Tx{N) = W(AA) = axN, (Al) 

where N is the column density of the absorbing species, / is the oscillator strength, (f) is the line profile 
function, and A A = A — Ao, where Ao is the wavelength at line center. The continuum reconstruction method 
finds the correct column density of an ionic species by multiplying the observed flux, Fx^ by exp[-|-TA(iV)]. 
The column density N is chosen such that the reconstructed continuum best matches the non-absorbed 
stellar continuum, Fx^o- 

This method is best applied to very strong lines with cores near zero flux (very large optical depths) 
and strong Lorcntzian damping wings. We rely on these wings to provide us with an accurate estimate of 
Tx{N). The Voigt line profile function is a convolution of natural and Doppler broadening components: 

0y ( AA) = ( AA) ® ( A A) . ( A2) 

The Lorentzian profile, (/)l(AA), is 

,^ ^, SA/47r2 ^ 

where is the summed Einstein probabilities for spontaneous radiative decay. Using the /-values and 
transition probabilities from Morton (1991), we derive the following numerical values in ax for the Lya 
absorption: 

, 2n f 4.26x10-20 1 , 

where A A is in A. The Dopper profile, 0£)(AA), is written 

0Z3(AA) - exp[-(cAA/A„6)2]. (A5) 

Here the Dopper spread parameter, b, includes thermal and non-thermal components added in quadrature 



[e.g., equation (10)]. To derive the thermal component we assume T — 6,000 K, equivalent to the 
temperature for component 1 derived in § [l.4| . For the non-thermal component we estimate the &- value of a 
single gaussian component fit to the S II Na(v) profile. Assuming the spread is primarily due to non-thermal 
broadenings, we derive u„j w 8 km s"^.^ The total Dopper spread parameter is 6 = 12.8 km s^^ [in good 
agreement with the value b = 12^2 derived by York & Rogerson (1976)]. 



^"For gas with a temperature of 6000 K, we expect a b-value b <2 km s ^ for S II, suggesting that the non-thermal motions 
do dominate the observed Na{v) profile of this ion. 
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To derive the reconstructed continuum flux, F^, we therefore perform the following 

pR _ ^ ( Ap,] 

^ exp[-aAiV]0 0/(AA)' ^ ' 

The additional term (pi {AX) in the denominator is the instrumental spread function. 

We adopt an instrumental spread function that consist of two Gaussian components: a narrow "core" 
and more extended "halo." The width of the core component was taken from the plots shown in the 
GHRS Instrument Handbook (Soderblom et al. 1995). The adopted resolutions of the G140M and G160M 
gratings at A1215.67 are Av = 14.4 km s~^ and 19.3 km s"^ (FWHM), respectively. The halo component 
of the adopted spread function is due to effects at the detector (see Cardelli, Ebbets, & Savage 1990). 
We characterise it as a Gaussian with a FHWM of 5 diodes (Spitzer & Fitzpatrick 1993). The relative 
contribution of these two components was derived following Spitzer & Fitzpatrick (1993) by fitting a spline 
to the relative contributions as a function of wavelength observed in pre-flight testing of the instrument 
(Cardelli et al. 1990). At A1215.67 the contribution of the halo component is 7.5% of the total spread 
function. 

The GHRS dataset for /z Col includes observations of the region containing interstellar Lya absorption 
(Ao = 1215.670) made through the SSA and LSA with both the G140M and G160M intermediate resolution 
gratings. To derive the H I column density we have used the G160M SSA data {S/N ^ 100) with the LSA 
and SSA G140M data {S/N ^ 60 and 30, respectively). We use observations taken through both apertures 
with the G140M data because the signal to noise of the SSA data do not allow us to probe as far into 
the core as we would like, while the LSA data are complicated by the poorly known instrumental spread 
function. The scattered light properties of the intermediate resolution GHRS first-order gratings are superb. 
Figure ^ shows the wavelength region near Lya observed with the G160M grating. The core of the deep 
Lya profile is fiat bottomed at zero flux level, consistent with no contribution from scattered light. The use 
of observations taken with both gratings reassures us that the complicating influences of spread functions, 
variations in the sensitivity at different points on the blaze function, and other instrumental effects do not 
impact our derived column density. In performing our analysis we follow the reconstructed continuum to a 
point where the observed signal to noise ratio falls below 10:1. Thus for the G160M observations we exclude 
information between 1215.0 < A < 1216.3 A. This corresponds to an optical depth (before instrumental 
smearing) of t ^ 5.0. The two observations with the G140M grating do not allow us to trace the continuum 
to quite as high optical depths. 

From our analysis of the GHRS observations of Lya towards fi Col, we derive a column density 
of neutral hydrogen logiV(H I) = 19.87 ± 0.015 (ilcr systematic). Figure |l^ shows the observed and 
reconstructed intensities for the G160M observations. We also show the ±2a (systematic) reconstructed 
continua. The derivation of A'^(H I) is relatively subjective. Given the complex undulating stellar continuum 
against which we are viewing the Lya absorption, we feel a formal fitting procedure is unwarranted. 

Estimating the errors in our column density determination is not a straightforward or easily 
quantifiable task. The undulations in the reconstructed continuum are primarily due to stellar absorption 
lines. The signal to noise ratios of our Lya spectra are high enough that the systematic errors of the 
continuum reconstruction method (i.e., continuum placement, central velocity, etc.) almost certainly 
dominate the statistical errors (e.g., photon statistics, background subtraction, etc.) in our derivation. 
We have estimated ±2a systematic errors for the values of A^(H I) we judged to give the best continuum 
reconstruction by increasing and decreasing A^(H I) until we judge the resulting continuum reconstructions 
as highly implausable. In Figure ^ we display the best fit continuum reconstruction along with the ±2ct 
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reconstructions. One can see that these ±2f7 results are truly unacceptable. That the range of appropriate 
values is so small is a result of the high signal to noise and exceptional scattered light properties of the 
observations. 

One problem in deriving the H I column density is determining the velocity alignment of the spectrum. 
We have used the lines of the N I triplet near A1200 (e.g., see Figure ||) to give us a range of acceptable 
values for the velocity zero-point of our spectrum, but we have used the reconstructed wings of the H I Lya 
profile themselves to fine-tune this determination. The assumption in this approach is that the absorption 
is symmetric. This may not be the case given complex component structure along the line of sight. Though 
the S II profile is reasonably approximated by a Gaussian fit (i.e., our turbulent broadening discussed 
above), we have in deriving our final number derived the thermal Gaussian profile with a detailed model of 
the component structure along the sightline. We use the results of our component fitting to the S II profile, 
with additional approximations for components 3, 4, and 5, as the kernel for the non-thermal structure 
along the sightline. The results of using this detailed model are almost indistinguishable from the simple 
Gaussian turbulent broadening model described above. 

Because the derivation of iV(H I) relies heavily upon the very broad wings of the absorption profile, 
the precise details of the convolutions described here make little difference to the final result. As mentioned 
above, using a detailed component model as opposed to a simple Gaussian model for the non-thermal 
broadening gives indistinguishable results. Adopting temperatures in the range 100 < T < 10000 K 
makes virtually no difference to the final continuum reconstruction. Even neglecting the smearing by the 
instrumental spread function makes only very negligible changes to the reconstructed continuum. Not 
only are the results we present quite robust to errors in these facets of our approach, but future works at 
comparable resolutions need not worry so much about the details of such complications when attempting 
to derive the H I column density along a given line of sight. Again, the keys to our ability to precisely 
derive the H I column density are the high signal-to-noise ratios of our data, the exceptional scattered light 
properties of the GHRS first-order gratings, and an appropriate (though not overly high) resolution. 

Our derived value of N(Il I) is quite similar to past determinations. From lUE measurements of Lya, 
Diplas & Savage (1994) derive logiV(H I) = 19.88±0.08. Bohlin et al. (1978) derive logiV(H I) = 19.85±0.08 
from Copernicus observations of Lya, and from Copernicus measurements of Lyman-/? and higher Lyman 
series absorption, York & Rogerson (1976) have estimated log-/V(H I) = 19.9 ± 0.1. The primary advantage 
of our derivation is that the errors are almost completely due to the systematics of the continuum 
reconstruction method applied to a star with undulating stellar absorption features ^ 5 — 30% deep. 

The derived column density logiV(HI) = 19.87 ± 0.015 contains contributions from the stellar 
photosphere as well as the interstellar medium. We must account for the stellar component if we are to 
very accurately derive the interstellar H I column density. Savage & Panek (1974) have shown that the 
observed equivalent widths of stellar Lya absorption for nearby early-type stars are in agreement with the 
non-LTE predictions of Mihalas (1972a, b). Their work was used by Diplas & Savage (1994) to derive stellar 
contributions to the Lya equivalent width based upon an empirical calibration using the Stromgren [ci] 
index. They suggest a contribution from the photosphere of /i Col of order ^^^(Lya) = 2.95 A, equivalent 
to logiVm = 18.27 + log M^A (Lya) = 18.74. 

In this work we will estimate the contribution of stellar Lya absorption by using the Mihalas (1972a, b) 
models presented in Savage & Panek (1974) to constrain the expected ratio of WxiLya) /Wx(ii0). Using 
measurements of the H/3 equivalent width by Buscombe (1969) we will then estimate the stellar contribution 
to WA(Lya). Using the data in Table 3 of Savage & Panek (1974) we find the ratio of these two hydrogen 
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lines is well fit by a third order polynomial such that: 

W,iLya)/W,m = 14.89- 38.33 (^5^) + 39.53 (5^)' - 15.12 {^^^ ■ (A7) 

For /i Col we predict W\{Lya) /W\{}iP) « 0.418. The measured 11/3 equivalent width from the photosphere 
of n Col is WxiUP) = 1.77 A (Buscombe 1969), suggesting a stellar contribution to the Lya equivalent 
width of Wx(Lya) ~ 0.74 A. This is equivalent to logiVni ~ 18.14. Thus the stellar contribution should 
represent « 0.01 dex of the total. We adopt a final interstellar column density to the star /i Col of 
logA^Hi = 19.86 ± 0.015 {la systematic). 

B. DERIVATION OF THE LARGE SCIENCE APERTURE LINE SPREAD FUNCTION 

In the process of attempting to fit component models to our interstellar line profiles, we have come 
to the conclusion that the LSF presented by Robinson et al. (1998) for the LSA is not appropriate at 
most wavelengths covered by the GHRS, evidenced by our inability to adequately fit model profiles to 
the narrow, often deep, component 4. Figure |l^ shows the best fit model for Mg II AA2796 and 2803 in 
the top row when using the LSF suggested by Robinson et al. (1998). Our best fit component models 
for Mg II at AA2796 and 2803 using the Robinson et al. LSF were unable to match the depth of the core 
of component 4, while at the same time they provided too much absorption in the wings of the profile. 
Furthermore, when we fit the pre-COSTAR observations of Mg II A2803 taken through the SSA using the 
spread function derived by Spitzer & Fitzpatrick (1993) and Fitzpatrick (priv. comm.), we found significant 
discrepancies in the derived column densities for components 3 and 4, which are both very narrow and 
deep. Using the Robinson et al. LSF we derived a column density logNiiMg II) = 13.22 ± 0.09; using the 
SSA observations and the LSF of Fitzpatrick & Spitzer we derived logiV4(Mg II) — 12.86 ± 0.05, more than 
a factor of two different. Also, the derived &- values were slightly different: 64 = 2.0 ± 0.3 km s^^ using the 
SSA observations and 64 = 1.54 ± 0.15 km s""'^ using the LSA observations. We found similar discrepancies 
comparing the SSA and LSA observations of O I. However, the SSA and LSA results agreed relatively well 
when comparing fits of the Si II A1304 transition. We believe this is because component 4 is not as deep in 
this line as the others. Thus the differences in the LSF seem to give significantly different results only for 
strong, very narrow absorption features. 

We have attempted to rederive the line spread function for the LSA. Though the problem is poorly 
constrained by the current dataset, we have pieced together a spread function that is able to match the 
SSA results and provide good fits to the data. To begin we have fit the Robinson et al. (1998) spread 
function with two Gaussians: a narrow "core" and a broad "halo" component. Table ^ gives the FWHMs 
(in diodes) and relative areas of these Gaussians. Our approach has been to assume a two-component 
Gaussian model is a valid approximation to the LSA spread function, and that the values of the FWHM 
from the Robinson et al. LSF are reasonable. We assume the widths of the core and halo Gaussians in 
diode-space are constant. The effective resolution changes with placement of the diode array in the echelle 
due to differences in the mapping of wavelengths to diode space. Our free parameter in the fit is the relative 
area of the core and halo components. We have proceeded by fitting the SSA data for O I, Si II, Mg II, and 
Fe II using the Fitzpatrick & Spitzer LSF. Then we have held the results of these fits constant and applied 
them to our LSA observations using varying values of the relative power of the core and halo components of 
our LSF. For each transition of these species that showed component 4 strongly, we have derived a best fit 
value for the relative power of the halo to the LSF. Figure |l^ shows our derived spread function with the 
Robinson et al. (1998) and Spitzer & Fitzpatrick (1993) LSFs at A1900 A. 
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The extended wings of the LSF in the LSA are a resuh of the broad wings of the point spread 
function from the aberrated primary (plus COSTAR) being aUowed to enter the spectrograph through the 
1'.'74 X 1'.'74 LSA. For the SSA, which effectively blocks the far wings of the telescope point spread function, 
the broad component adopted by Fitzpatrick & Spitzer is a result of aberrations in the internal optics of 
the spectrograph and effects in the detector. The relative power of the broad component is much greater 
in the spread function for the LSA than for the SSA. We have fit the fractional contribution of the halo 
component, fhaio, to the total as a function of wavelength with the following function: 

= 0.486 - (2.89 X 10"'*)A + (5.52 x 10^*)A^ (Bl) 

We find that this prescription gives results that agree in column density and 5-value with fits to the SSA 
data using the Fitzpatrick & Spitzer SSA spread function. The best fit profile using this profile is also better 
able to reproduce the data than the Robinson et al. LSF. Figure |8| shows the results of the component fits 
for the Mg II lines. The top row shows the best fit component model using the Robinson spread function 
for the Mg II lines near A2800, the middle row shows the best fit using our derived spread function, and the 
bottom row shows the fit to the SSA data using the Spitzer & Fitzpatrick (1993) spread function. Also 
shown in the bottom row is the fit obtained simultaneously for the Mg II A1239 transition using our derived 
LSF. The tick marks show the locations of our model interstellar components. When examining the fit to 
component 4 near wlsr = +41 km s~^, it is clear that our LSF provides a better fit to this strong, narrow 
component. Furthermore, it provides results consistent with those derived using the SSA data. 

While our empirically-derived LSF is able to match the SSA data, we do not wish to suggest the LSA 
spread function is well understood. Our fit is consistent with the available data, though the problem is 
ill-constrained, particularly in the wavelength range 1600 < A < 2300 A, where we have no data. The true 
LSF of the LSA is likely more complex than the two Gaussian model we have adopted. Further, the true 
spread function may be a function of the positioning of the star within the LSA. It is important to realize 
that the adoption of a specific shape for the instrumental profile can have important consequences when 
trying to derive accurate column densities for narrow, deep lines. When possible, investigators interested 
in deriving the column densities of such components should use the SSA, for which the spread function 
has been better characterized (e.g., Spitzer & Fitzpatrick 1993) and is better behaved. When using data 
taken through the LSA for such work, one should investigate the consequences of adopting slightly different 
spread functions. In our case the values derived for the column densities of O I and Mg II were different by 
factors of 1.5 to 2.5 between the different spread functions. 
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Table 1. Velocities and Conditions in the Absorbing Regions Towards fj, Columbae 



Region 






n~i c 

e 




Previous 


Number 


[km s-i] 


[km s""'^] 


[K] 


[cm~^] 


ID 


1 


3 


-17.0 to 15.5 


6, 200 ±300 


0.27 ±0.05 


A 


2 


20.1 


15.5 to 29.0 


6, 900 ± 600 


0.31 ±0.13 


B 


3 


31.0 


29.0 to 37.5 


^10,000 




D 


4 


41.2 


37.5 to 47.0 


4, 000 ± 700 


0.64 ±0.12 


C 


5 


-28.8 


-45.0 to -17.0 


^ 36,000 




... 6 



Approximate central velocity of the absorbing region in the LSR frame. For 

components 2, 3, and 4 these values come from fitting a model interstellar line 
profile to the Si II lines, which are absolutely calibrated using the nearby O I* 
atmospheric line. 

''The velocity range over which we have integrated to obtain the gross values of 
equivalent width and apparent column density for these absorbing regions. 

'^Electron temperatures as derived from ionization equilibrium and component 
fitting analyses. For components 1 and 2, these values were derived from ionization 
equilibrium of Ca and Mg (§4.4.2). For components 3, 4 and 5 the temperatures 
are derived from our component fitting results (§5.1). 

•^Electron densities as derived from various combinations of ionization 
equilibrium. For components 1 and 2, these values were derived from ionization 
equilibrium of Ca and Mg (§4.4.2). For component 4 these values were derived from 
ionization equilibrium of Mg assuming the temperature suggested from component 
fitting methods (§5.1). 

"ShuU & York (1977) refer to this component as "trailing absorption" in their 
discussion of the Si III profile. 



Table 2. Log of GHRS Observations 



Spectral Range 


Rootname** Exp.'' Mode Aper.^ dc** 


Species 


[A] 


[sec] & Order'^ 


Detected 



1116- 


■1124 


Z2X30108T 


81.6 


G140M/1 


LSA 


0.0 


Fe II, Fe III 






Z2X30109T 


81.6 


G140M/1 


LSA 






1131- 


■1137 


Z2C0011AP 


108.8 


A/49 


LSA 


0.011 


N I, Fe II 


1142- 


■1148 


Z2AF011AT 


108.8 


A/49 


LSA 


0.061: 


Fell 






Z2C0030RT 


108.8 


A/49 


LSA 










Z308010OT 


108.8 


A/49 


LSA 










Z3JN010OT 


108.8 


A/49 


LSA 






1149- 


■1153 


Z2AF020RT 


108.8 


A/49 


LSA 


0.061: 


PII 






Z2C0010QP 


108.8 


A/49 


LSA 










Z2C0021AP 


108.8 


A/49 


LSA 






1154- 


■1161 


Z2C00119P 


108.8 


A/48 


LSA 


0.058: 


CI 


1157- 


■1163 


Z2C0041AT 


108.8 


A/49 


LSA 


0.011: 


CI 


1184- 


■1191 


Z2AF010PT 


108.8 


A/47 


LSA 


0.000 


C I, Si II, S III 






Z2C0020PP 


108.8 


A/47 


LSA 






1191- 


■1197 


Z2AF0118T 


108.8 


A/47 


LSA 


0.038 


Si II 






Z2AF021CT 


108.8 


A/47 


SSA 










Z2C0030PT 


108.8 


A/47 


LSA 






1195- 


•1230 


Z2AN020AT 


54.4 


G160M/1 


SSA 


0.0 


HI 


1197- 


■1225 


Z2XE010DT 


27.2 


G140M/1 


LSA 


0.0 


HI 


1197- 


-1225 


Z2XE0117T 


27.2 


G140M/1 


SSA 


0.0 


HI 


1197- 


■1200 


Z2AF020PT 


108.8 


A/47 


LSA 


0.038: 


N I, Mn II 


1200- 


■1203 


Z2C00218P 


108.8 


A/47 


LSA 


0.035 


NI 


1206- 


■1212 


Z2AZ010MTf 


1740.8 


A/46 


SSA 


0.027 


Si III 


1237- 


■1244 


Z2AZ010CT^ 


1088.0 


A/45 


LSA 


0.049: 


N V, Mg II 






Z2AZ010OTf 


1088.0 


A/45 


SSA 






1249- 


■1256 


Z2AF020MT 


27.2 


A/45 


LSA 


0.049: 


SII 






Z2C00215P 


27.2 


A/45 


LSA 






1259- 


■1262 


Z2C0040MT 


27.2 


A/45 


LSA 


0.033 


SII 


1259- 


■1267 


Z2C00114P 


27.2 


A/44 


LSA 


0.026 


Si ir 


1271- 


■1278 


Z2AF0114T 


27.2 


A/44 


LSA 


0.046: 


CI 


1301- 


■1308 


Z2AF0113T 


27.2 


A/44 


LSA 


0.050 


I, Si II, P II 






Z2AF0217T 


27.2 


A/44 


LSA 










Z2C0030KT 


27.2 


A/44 


LSA 










Z308010HT 


27.2 


A/44 


LSA 






1313- 


■1320 


Z2C0010JP 


27.2 


A/43 


LSA 


0.050: 


Ni II 






Z2C0040KT 


27.2 


A/43 


LSA 






1326- 


■1333 


Z2AF010JT 


27.2 


A/42 


LSA 


0.041: 


CI 






Z2C0020JP 


27.2 


A/42 


LSA 






1332- 


■1339 


Z2AF0112T 


27.2 


A/42 


LSA 


0.041 


c II, c ir 






Z2AF0216T 


27.2 


A/42 


SSA 










Z2C0030JT 


27.2 


A/42 


LSA 
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Table 2 — Continued 



Spectral Range 


Rootname'* 


Exp.i^ 


Mode 


Aper."^ 


dc'^ 


Species 


[A] 




[sec] 


& Order^ 






Detected 




Z308010GT 


27.2 


A/42 


LSA 






1364-1372 


Z2AF0111T 


27.2 


A /At 

A/41 


LSA 


0.039 


Ni II 




Z2AF0215T 


27.2 


A/41 


SSA 








Z2C0030IT 


27.2 


A /At 

A/41 


LSA 








Z308010FT 


27.2 


A /At 

A/41 


T O A 

LSA 






1526-1534 


/2C0040DT 


tr A A 

54.4 


A/37 


T O A 

LSA 


0.055 


Si II 




Z2U00401j1 


27.2 


A/37 


T O A 

LSA 






1554-1562 


Z2AF010VT 


108.8 


A/36 


LSA 


0.031 


C I 




Z2AF020ZT 


108.8 


A/36 


SSA 








Z2C0030CT 


54.4 


A/36 


LSA 






1606-1614 


Z2AF020BT 


108.8 


A/35 


LSA 


0.030 


Fe II 




Z2C0020UP 


54.4 


A/35 


LSA 






1653-1662 


Z2AF020AT 


108.8 


A/34 


LSA 


0.030 


C I 




Z2C0020TP 


54.4 


A/34 


LSA 






1668-1676 


Z2C0040TT 


54.4 


A/34 


LSA 


0.030 


Al II 


1706-1712 


Z2AF0209T 


108.8 


A/33 


T n A 

LSA 


0.029 


Ni II 




Z2CX011KT 


27.2 


B/33 


T O A 

LSA 


0.045 






Z2C0020bP 


54.4 


A/33 


LSA 


0.029 






Z2D4021JT 


54.4 


B/33 


LSA 


0.045 




1740-1749 


Z2D4020MT 


54.4 


B/32 


T O A 

LSA 


0.043 


AT* TT 

Ni 11 


1 D n f 1 D (\C\ 

1805-1809 


Z2CX010Mi 


27.2 


B/31 


LSA 


0.041 


C TT 

Si II 




Z2D40211 1 


54.4 


B/31 


T O A 

LSA 






1846-1857 


Z2D40118T 


54.4 


B/30 


LSA 


0.039 


Al III 


1860-1866 


Z2Ca010LT 


27.2 


B/30 


LSA 


0.039 


Al III 




Z2D4020KT 


54.4 


B/30 


LSA 






2017-2027 


Z2D4030IT 


54.4 


B/28 


LSA 


0.035 


Mg I, Zn II 


2052-2061 


Z2D40115T 


54.4 


B/27 


LSA 


0.033 


Cr II 


2062-2063 


Z2D40115T 


54.4 


B/27 


LSA 


0.033 


TT 

Cr II 




Z2D4020HT 


54.4 


B/27 


LSA 






2063-2074 


Z2D4020HT 


54.4 


B/27 


LSA 


0.033 


Cr II 


2242-2254 


Z23L040FT 


57.6 


B/25 


LSA 


0.030 


Fe II 




ZiZyiUbUr i 


57.0 


rJ/zo 


T O A 

LbA 








Z2D4021BT 


54.4 


B/25 


LSA 








Z308020DT 


115.2 


B/25 


LSA 








Z3JN020DT 


115.2 


B/25 


LSA 






2259-2264 


Z2CX011CT 


27.2 


B/25 


LSA 


0.030 


Fe II 




Z2D4030FT 


54.4 


B/25 


LSA 






2335-2348 


Z23L040ET 


57.6 


B/24 


LSA 


0.031 


Fe II 




Z308020CT 


115.2 


B/24 


LSA 








Z3JN020CT 


115.2 


B/24 


LSA 






2364-2376 


Z2D4030UT 


54.4 


B/24 


LSA 


0.028 


Fe II 
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Table 2 — Continued 



Spectral Range 


Rootname^ 


Exp> 


Mode 


Aper.'^ 


do*" 


Species 


[M 

L J 




[seel 

L J 


& Order'^ 






Detected 


2567-2581 


Z2D4030CT 


54.4 


B/22 


LSA 


0.017: 


Mn II 


2575-2590 


Z2D4030ST 


54.4 


B/22 


LSA 


0.017: 


Fell 


2590-2603 


Z2D40318T 


54.4 


B/22 


LSA 


0.017 


Mn II, Fe II 


2792-2802 


Z23L0413T 


28.8 


B/20 


LSA 


0.015 


Mg II 




Z29I0613T 


28.8 


B/20 


LSA 






2802-2807 


Z23L0413T 


28.8 


B/20 


LSA 


0.015 


Mg II 




Z23L0414T 


28.8 


B/20 


LSA 








Z23L040AT 


57.6 


B/20 


LSA 








Z29I060AT 


57.6 


B/20 


LSA 








Z29I0613T 


28.8 


B/20 


LSA 








Z29I0614T 


28.8 


B/20 


LSA 






2852-2857 


Z23L0418T 


28.8 


B/20 


LSA 


0.015: 


Mgl 




Z23L0419T 


28.8 


B/20 


LSA 








Z29I0618T 


28.8 


B/20 


LSA 







'^STScI archival rootname. 

''Total exposure time given in seconds. Typically 94% of this time was spent on- 
spectrum, while 6% was used to measure the interorder echelle background. 

'^Echelle mode and spectral order used for the observation. 

^Aperture used for the observation. The LSA subtends 1'.'74 x 1'.'76 on the sky, the 
SSA 0'.'22 X 0'.'22, for post-COSTAR observations. 

"^Second-order scattered light correction (see Cardelli ct al. 1993). For spectral 
regions containing strongly saturated lines we have derived dc directly. Those values 
marked with a colon are regions for which we have not independently confirmed the 
value of dc- In several cases we have applied our derived values of dc to spectral 
regions in the same echelle order that lack strongly saturated lines; in regions where no 
independent information was available we have used the values suggested by Cardelli 
et al. (1993) for pre-COSTAR SSA observations. 

^This observation was taken using the FP-SPLIT=4 observing routine. The fixed- 
pattern noise spectrum of this observation has been removed from the observations. 



Table 3. Integrated Equivalent Widths 



IP IP Wa [mA] 

Species x'"^ X' Kac"" log A/ Ref.'^ Comp 1 Comp 2 Comp 3 Comp 4 S/N^ 

[eV] [eV] [A] (3 km s'^) (20 km s"!) (31 km s'^) (41 km s'^) 



CI 




11.26 


1328.833 


1.887 


1 


7.2 ±0.5 


< 2.2« 


< 1.9 


< 1.9 


48 








1560.309 


2.099 


1 


12.0 ±0.7 


2.7 ±0.5 


< 1.7 


< 1.7 


64 








1656.928 


2.367 


1 


18.8 ±2.4 


3.6 ±1.5 


< 3 


< 3 


26 


cr 




11.26 


1656.267 


1.987 


1 


2.4 ±0.9 


< 3 


< 3 


< 3 


42 








1657.907 


1.890 


1 


4.8 ±2.0 


<4 


< 4 


<4 


32 


cr* 




11.26 


1277.550 


2.016 


1 


<3 


< 2.0 


< 1.7 


< 1.7 


51 








1561.438 


2.023 


1 


< 4 


< 3 


< 3 


< 3 


42 








1658.121 


1.765 


1 


< 5 


<4 


< 4 


<4 


31 


CII 


11.26 


24.38 


1334.532 


2.232 


1 


155 ±9 


56.0 ± 1.1 


65.0 ±1.3 


28.3 ±0.6 


56 


Cn* 


11.26 


24.38 


1335.708 


2.186 


1 


132 ±3 


23.6 ±0.6 


2.3 ±0.6 


3.5 ±0.5 


31 


NI 




14.53 


1134.165 


1.182 


1 


79.2 ±2.0 


28.0 ± 1.0 


< 3 


< 3 


12 








1134.415 


1.483 


1 


96.2 ±2.1 


24.9 ± 1.0 


< 3.3 


2.9 ± 1.0 


12 








1134.980 


1.660 


1 


106.0 ±2.2 


23.8 ±0.9 


< 2.4 


5.0 ±0.8 


14 








1199.550 


2.202 


1 


118 ±3 


34.9 ±0.8 


f 


14.8 ±0.8 


27 








1200.223 


2.026 


1 


116.0 ±2.3 


34.3 ±0.7 


f 


11.3 ±0.4 


46 








1200.710 


1.725 


1 


109.2 ±2.2 


32.2 ±0.7 


f 


7.7 ± 0.4*^ 


49 


01 




13.62 


1302.169 


1.804 


1 


137 ±3 


52.6 ± 1.1 


16.4 ±0.4 


22.8 ±0.5 


54 








1355.598 


-2.772 


1 


< 3 


< 2.3 


< 2.0 


< 2.0 


46 


Mgl 




7.65 


2026.477 


2.3698 


1 


16.5 ± 1.7 


2.0 ±0.5 


< 1.5 


< 1.5 


44 








2852.964 


3.718S 


1 


207 ±4 


44.8 ± 1.0 


2.2 ±0.4 


3.9 ±0.5 


80 


Mg II 


7.65 


15.04 


1239.925 


-O.lOOe 


2 


9.82 ±0.24 


0.79 ± 0.09 


< 0.3 


< 0.3 


195 








1240.395 


-0.401S 


2 


5.47 ±0.18 


0.36 ±0.09 


< 0.18 


< 0.18 


217 








2796.352 


3.234 


1 


307 ±6 


117 ±3 


59.4 ± 1.2 


49.9 ± 1.1 


111 








2803.531 


2.933 


1 


297 ±6 


112.1 ±2.3 


46.3 ±1.0 


44.2 ± 0.9 


136 


Al II 


5.99 


18.83 


1670.787 


3.486 


1 


148 ±3 


55.7 ± 1.3 


9.3 ±1.0 


6.5 ± 1.0 


16 


Al III 


18.83 


28.45 


1854.716 


3.017 


1 


15.0 ±2.0 


< 3.0 


< 3.0 


< 3.0 


30 








1862.790 


2.716 


1 


9.9 ± 1.5 


< 2.4 


< 2.4 


< 2.4 


35 



Table 3 — Continued 



IP IP Wa [mA] 

Species x'"^ A^ac*" log A/ Ref> Comp 1 Comp 2 Comp 3 Comp 4 S/N^^ 

[eV] [eV] [A] (3 km s'^) (20 km s'^) (31 km s'^) (41 km s'^) 



Si 11 


8.15 


16.35 


1190.416 


2.474 


1 


116.0 ±2.3 


48.3 ± 1.0 


12.1 ± 0.3 


12.3 ±0.3 


54 








1193.290 


2.775 


1 


124 ±3 


50.0 ± 1.0 


19.0 ±0.5 


17.0 ±0.4 


48 








1304.370 


2.050 


3,4 


121 ±3 


42.5 ± 0.9 


6.91 ± 0.24 


8.5 ± 0.3 


57 








1526.707 


2.225 


3,4 


148 ±3 


56.3 ± 1.2 


11.9 ±0.5 


12.4 ±0.6 


27 








1808.013 


0.596 


5 4 


61.2 ± 1.6 


6.5 ± 0.2 


< 1.5 


< 1.5 


36 


Si 11* 


8.15 


16.35 


1264.738 


3.058 


1 


3.60 ± 1.16 


< 1.8 


< 1.8 


< 1.8 


33 


Si 111 


16.35 


33.49 


1206.500 


3.304 


1 


126 ±3 


40.9 ± 0.9 


10.7 ±0.5^ 


1.4 ±0.5^ 


76 


P 11 


10.49 


19.72 


1152.818 


2.435*5 


1 


53.2 ± 1.4 


4.9 ±0.5 


< 1.5 


< 1.5 


32 








1301.874 


1.351*5 


1 


6.6 ± 0.6 


h 


h 


h 


55 


S 11 


10.36 


23.33 


1250.534 


0.834 


1 


68.6 ± 1.4 


7.4 ± 0.4 


< 0.9 


< 0.9 


40 








1253.811 


1.135 


1 


90.9 ± 1.8 


11.3 ±0.4 


< 0.9 


< 0.9 


43 








1259.519 


1.311 


1 


104.6 ±2.1 


20.0 ± 0.6 


< 1.3 


< 1.2 


36 


S 111 


23.33 


34.83 


1190.208 


1.421 


1 


16.3 ± 0.6 


< 0.9 


< 0.9 


h 


58 


Cr 11 


6.77 


16.50 


2056.254 


2.334 


6 


12.2 ± 1.3 


4.2 ±0.7 


< 1.5 


< 1.5 


43 








2062.234 


2.206 


6 


6.7 ± 1.3 


3.5 ± 0.6 


< 1.2 


< 1.2 


73 








2066.161 


2.027 


6 


6.0 ± 0.9 


2.6 ± 0.5 


< 1.2 


< 1.2 


48 


A/Tn TT 

iVlli ii 






1 1 07 1 84 




1 

L 




1 n -1- n 

1 .u HI u.o 


\ u.u 


\ u.u 


49 








2576.877 


2.956 


1 


55.8 ±1.3 


10.9 ±0.4 


1.6 ±0.4 


1.1 ±0.4 


68 








2594.499 


2.847 


1 


43.6 ± 1.1 


8.7 ±0.5 


< 1.2 


< 1.2 


60 


Fell 


7.87 


16.16 


1133.665 


-1.2228 


1 


9.5 ±2.2 


4.5 ±1.1 


< 3 


< 3 


13 








1143.226 


1.1828 


1 


24.7 ± 1.0 


10.1 ±0.4 


< 2.2 


< 2.2 


35 








1144.938 


2.080 


1 


75.1 ±1.7 


32.1 ±0.8 


3.6 ±0.6 


2.1 ±0.6 


31 








1608.451 


1.998 


7,8 


100.3 ±2.1 


36.6 ±0.8 


3.6 ±0.5 


3.9 ±0.5 


35 








2249.877 


0.612 


9,10 


10.2 ±0.3 


3.96 ±0.20 


< 0.6 


< 0.6 


218 








2260.780 


0.742 


9,10 


15.0 ±0.5 


6.9 ±0.3 


< 0.9 


< 0.9 


97 








2344.214 


2.4108 


1 


204 ±4 


70.0 ± 1.4 


16.0 ±0.3 


11.7 ±0.3 


236 








2374.461 


1.8898 


7,10 


125 ±3 


47.7 ± 1.0 


5.3 ±0.3 


4.1 ±0.3 


65 








2586.650 


2.248 


7,10 


196 ±4 


74.0 ± 1.5 


11.4±0.4 


9.8 ±0.4 


60 








2600.173 


2.765 


1 


238 ±5 


98.9 ±2.0 


28.7 ±0.7 


22.1 ±0.6 


51 



Table 3 — Continued 





IP 


IP 










Wa 


[mA] <= 






Species 




X' 


\ ^ 


log A/ 


Ref.b 


Comp 1 


Comp 2 


Comp 3 


Comp 4 


S/Nd 




[eV] 


[eV] 


[A] 






(3 km s-i) 


(20 km s-i) 


(31 km s-i) 


(41 km s-i) 




Felll 


16.16 


30.65 


1122.526 


1.947 


1 


19 ±4' 








24 


m; tt 


1.04: 


10 17 


1 Q1 7 01 7 

ioi / .Zi 




111 

1,11 


c -1- n 7 
D.S ± U. / 


6.0 ± U.4 


< i.z 


< i.z 


A A 

44 








1370.132 


2.023s 


11,12 


8.1 ±0.5 


3.4 ± 0.3 


< 0.6 


< 0.6 


68 








1709.600 


1.784 


12 


4.9 ±1.0 


3.4 ±0.5 


< 1.2 


< 1.2 


44 








1741.569 


1.858 


12 


14.5 ±1.6 


5.5 ±0.9 


< 2.4 


< 2.4 


21 


Cu II 


7.73 


20.29 


1358.773 


2.713 


1 


< 3 


< 2.2 


< 1.9 


< 1.9 


50 


Zn II 


9.39 


17.96 


2026.136 


2.9966 


6 


52.5 ±1.6 


4.5 ±0.8 


... g 


... g 


46 








20(i2.G(i4 


2.72;-!K 


(i 


2(i..', ± 1.1 


2.2 ± ()..") 


< 1.2 


< 1.2 


Go 



'^Vacuum wavelengths from Morton (1991). 
'^Reference codes for /—values. 

■^Measured equivalent widths and ±lo- errors, given in mA. 
^Empirical signal-to-noise ratio derived from the fit to the continuum. 
''AH upper limits are 3 a limits. 

^The values derived for these components are rather uncertain due to the nature of the continuum placement and the 

possibility of fixed-pattern noise contamination (see text). 

^Denotes oscillator strengths that wc feel may have uncertainties approaching or exceeding our measurement errors. This 
classification is based on the current dataset and our experience with the current oscillator strength measurements and other 
absorption line datasets. 

''Confused with adjacent line. 

'Based upon G140M data. This line has been integrated over the velocity range -20 km s~^ to +23 km s~^. Though the 
line C I Al 122.518 A potentially overlaps the Fc III absorption, wc believe its contribution to the value of Wx quoted here is 
negligible based upon the absence of the stronger line of C I at 1122.438 A. 

References. — (1) Morton (1991); (2) Fitzpatrick (1997); (3) Spitzer & Fitzpatrick (1993); (4) Dufton, et al. (1983, 1992); 
(5) Bergeson & Lawlcr (1993a); (6) Bergeson & Lawlcr (1993b); (7) CardeUi & Savage (1995); (8) Bergeson, et al. (1996a); 
(9) Bergeson, et al. (1994); (10) Bergeson, et al. (1996b); (11) Fedchak & Lawler (1999); (12) Zsargo & Federman (1997) 
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Table 4. Integrated Apparent Column Densities 



CII 
CII* 



Log(Na) [cm 



-21 b 



Species 


\ a 

[A] 


Comp 1 
(3 km s-i) 


Comp 2 
(20 km s-i) 


Comp 3 
(31 km s-^} 


Comp 4 
(41 km s-i) 


CI 
CI 
CI 

C I* 
CI* 


1328.833 
1560.309 
1656.928 
1656.267 
1657.907 


12.92 ±0.04 
12.87 ±0.03 
12.79 ±0.07 

12.651°:^! 


12.21 ±0.09 







1334.532 
1335.708 



X X x"^ 

XXX 



XXX 

> 13.3 



XXX 

12.13 ±0.10 



> 13.5 
12.32 ±0.14 



NI 


1134.165 


X 


X 


X 


X 


X 


X 




NI 


1134.415 


X 


X 


X 


X 


X 


X 


13.02^0:13 


NI 


1134.980 


X 


X 


X 


X 


X 


X 


13.08 ±0.08 


NI 


1199.550 


X 


X 


X 


X 


X 


X 


12.95 ±0.10^^ 


NI 


1200.223 


X 


X 


X 


X 


X 


X 


13.09 ±0.02 


NI 


1200.710 


X 


X 


X 


X 


X 


X 


13.19 ±0.10<^ 



01 



1302.169 



XXX 



XXX 



13.49 ±0.02 13.72 ±0.02 



Mgl 2026.477 12.63 ± 0.05 11.72 ±0.11 
Mgl 2852.964 12.54 ±0.02 11.68 ±0.02 



Mg II 


1239 


925 


Mg II 


1240 


395 


Mg II 


2796 


352 


Mg II 


2803 


531 


Al II 


1670 


787 


Al III 


1854 


716 


Al III 


1862 


790 


Si II 


1190 


416 


Si II 


1193 


290 


Si II 


1304 


370 


Si II 


1526 


707 


Si II 


1808 


013 


Si ir 


1264 


738 


Si III 


1206 


500 



p II 



15.08 ±0.02 

15.11 ±0.03 

XXX 
XXX 

XXX 

11.97 ±0.07 
12.08 ±0.07 

XXX 
XXX 

XXX 
XXX 

15.10 ±0.02 
11 47+0.12 



1152.818 13.46 ±0.02 



13.96 ±0.05 



13.92 



-0.10 



-0.12 
XXX 

XXX 



XXX 



14.05 ±0.06 

> 12.67 
12.29 ±0.05 



10.23 ±0.09 10.48 ±0.06 



12.44 ±0.02 12.35 ±0.02 
12.52 ±0.02 12.53 ±0.02 

11.30 ±0.05 11.21 ±0.06 



X 


X 


X 


12.70 ±0.02 


12.70 ±0.03 


X 


X 


X 


12.70 ±0.02 


12.60 ±0.02 


X 


X 


X 


12.79 ±0.03 


12.87 ±0.03 


X 


X 


X 


12.80 ±0.03 


12.82 ±0.03 



11.81 ±0.03*= 
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Table 4 — Continued 



Log(N„) [cm 



-21 b 



Species 


\ ^ 
^vac 


Comp 1 


Comp 2 


Comp 3 


Comp 4 




[A] 


(3 km s-i) 


(20 km s-i) 


(31 km s-i) 


(41 km s-i) 


PII 


1301.874 


13.42 ±0.04 








SII 


1250.534 


15.19 ±0.02 


14.04 ±0.03 






SII 


1253.811 


15.17 ±0.02 


13.97 ±0.03 






SII 


1259.519 


15.13 ±0.02 


14.08 ±0.03 







SIII 1190.208 13.82 ±0.02 



Cr II 


2056, 


.254 


Cr II 


2062, 


.234 


Cr II 


2066, 


.161 


Mn II 


1197, 


.164 


Mn II 


2576, 


.877 


Mn II 


2594, 


.499 


Fell 


1133, 


.665 


Fell 


1143, 


.226 


Fell 


1144, 


.938 


Fell 


1608, 


.451 


Fell 


2249, 


.877 


Fc II 


2260, 


.780 


Fell 


2344, 


,214 


Fell 


2374, 


.461 


Fc II 


2586, 


.650 


Fell 


2600, 


.173 


Felll 


1122, 


.526 


Ni II 


1317, 


.217 


Ni II 


1370, 


.132 


Ni II 


1709, 


.600 


Ni II 


1741, 


.569 


Zn II 


2026 


.136 


Zn LI 


20(i2 


,(i()4 



12.51 ±0.05 12.05 ±0.07 

12.37 ±0.09 

12.50 ±0.07 12.141^:^^ 



12.44 ±0.04 
12.50 ±0.02 
12.48 ±0.02 11.75 ±0.03 



12.50 ± 0.02 11.75 ±0.03 10.89^^2 WmtoH 



14.18 ±0.09 13.85^o:l2 

14.27 ±0.07 13.87 ±0.05 

14.15 ±0.03 13.70 ±0.02 

14.10 ±0.02 13.63 ±0.02 12.43 ±0.06 12.46 ±0.07 

14.11 ±0.02 13.70 ±0.02 

14.15 ±0.02 13.81 ±0.02 

XXX 13.52 ±0.02 12.55 ±0.02 12.41 ± 0.02 

14.09 ±0.03 13.67 ±0.02 12.54 ±0.04 12.43 ±0.04 

14.10 ±0.02 13.61 ±0.02 12.51 ± 0.03 12.43 ±0.03 
XXX XXX 12.50 ±0.02 12.35 ±0.02 



13.37tnf 



12.77 ±0.05 12.50 ±0.05 

12.82 ±0.03 12.46 ±0.04 

12.73 ±0.09 12.58 ±0.07 

13.16 ±0.06 12.71 ±0.07 

12.56 ±0.03 11.35 ±0.09 

12.4<S±0.0.S ll..S8±0.n9 



''Vacuum wavelengths from Morton (1991). 

^Integrated column densities for the components considered here with ±1(7 errors, 
given in atoms cm^^. 

'^The XXX indicates the profile is strongly saturated in this wavelength region, 
yielding no useful information. 

•^The Al 199.550 A and A1200.710 A lines of N I each show evidence for possible 

contamination from fixed-pattern noise as well as lying in regions where the continuum 
placement is uncertain. The errors quoted here are estimates based upon experiments 
with various reasonable continuum fits. 



Table 5. Component Fitting Column Densities and 6- values 





Comp 1 


Comp 2 


Comp 3 


Comp 4 


Comp 5 


Species 
X 


\ogN{XY 
[cm~^] 


¥- 

[km s~^] 


logiV(X)'* 
[cm~^] 


[km s~-^] 


logA^(X)'^ 
[cm~2] 


[km s-i] 


log TV (X)'* 
[cm~^] 


b" 

[km s~-^] 


log Ni^X)"^ 
[cm~^] 


ClI 


XXX 


XXX 


XXX 


XXX 


XXX 


o.U It U.O 


1 Q 7c; -L n nt^ 

lo. / It U.UO 


7 n -1- n Q 
1 .u It u.y 


1 Q Q7 -L n 1 n 
lo.o 1 It u.iu 


CII* 


X X x** 


o.O It U.O 


lo.ZO It U.Uo 






9 4 _i_ 1 c 

It 1.0 


IZ.oO It U.UO 






NI 


XXX 


A ^ -\- (\ P, 

fi. i It U.O 


1 A or _i_ n no 
±4t.Z0 It U.Uo 






9 1 7 -L n 9n 
Z.l / It U.ZU 


1 Q 1 7 -L n 09 
lo.l 1 It u.uz 






01 


XXX 


o.y It U.O 


14.00 It U.UO 


4 9-1-07 
4.Z It U. ( 


lo.Oo It U.Uo 


9 t^S -1- n 1 9 
Z.0(5 It U. IZ 


1 Q oq I n ri9 
lo.oo It U.UZ 






Mgl 


12.59 ±0.05 


Z.DU It u. lo 


11 79 -1- n 09 
L L. 1 Z It U.UZ 




1 n 98 -t- n 11 

lU.Zo It U. 1 1 




1 n t;9 -L n nzt 
lU.OZ It U.U4 






Mg II 


15.09 ±0.20 


/I Q _i_ 1 c; 

^.y It ±.0 


lo.yo It u.iz 


Z.Q It U.O 


19 en -1- n 09 
iZ.OU It u.uz 


1 89 -1- n no 
i.oz It u.uy 


1 9 S9 -L n 09 

iz.oz It u.uz 


7 -1- 1 1 
/ .0 It 1.1 


1 1 Q8 -1- n n9 
1 1 .yo It u.uz 


Al II 


XXX 


3.5 ±0.8 


12.64 ±0.05 


2.65^ 


11.47 ±0.05 


1.74^ 


11.27 ±0.06 






Sill 


15.13 ±0.11 


3.6 ±0.6 


13.98 ±0.05 


2.7 ±0.4 


12.78 ±0.02 


1.67 ±0.09 


12.85 ±0.02 




11.64 ±0.10 


Q ^ T TT 
Ol ill 


XXX 


4.6 ±2.0 


12.65 ±0.08 


3.34 ± 1.7 


11.83 ±0.10 




10.67 ±0.22 


7.5 ± 1.0 


12.29 ±0.03 


PII 


13.50 ±0.15 


3. 68 


11.98 ±0.13 














SII 


15.24 ± 0.05 


5.0 ±0.7 


13.82 ±0.02 












1 

QJl 


Cr II 


12.48 ±0.07 




12.03 ±0.05 












CO 

1 


Mn II 


12.52 ±0.10 


3.0 ± 1.0 


11.61 ±0.05 


1.76' 


10.82 ±0.09 


1.18' 


10.55 ±0.14 






Fell 


14.13 ±0.02 


2.75 ±0.07 


13.78 ± 0.02 


1.76 ±0.14 


12.63 ±0.02 


1.18 ±0.11 


12.61 ±0.02 


o q+50 


11.34 ±0.10 


Ni II 


12.86 ±0.08 


2.79j 


12.35 ±0.04 














Zn II 


12.55 ±0.08 


2.79j 


11.16 ±0.08 















''Total column density contained within the component 1 blend. We have fit this region with three individual components, but we choose only 
to display the sums given the lack of uniqueness often involved with fitting this central region. The values given here should be compared with the 
values of the integrated column densities given in Table 4. 

"^Dopplcr spread parameters and ztla errors (in km s~^) derived through component fitting for the absorbing components 2, 3, and 4. These values 
were derived through component fitting techniques as described in §3.3.2. 

■^Column densities and ±1(7 errors (in atoms cm^^) derived through component fitting for the absorbing components 2, 3, and 4. 

*The XXX indicates the profile is strongly saturated in this wavelength region, yielding little or no reliable information. 

*The 6-values adopted for components 3 and 4 in the Al II fit are the weighted means of the Mg II and Si II fits for the same components. 

^The 6- values adopted for component 2 in P II were set equal to that derived for Si II: 6 = 3.6 ± 0.6 km s~^. 

^The 6-values for component 2 in this ion were set equal to that derived for Mn II: b = 3.0 ± 1.0 km s~^. 

'The lvalues for components 3 and 4 for the Mn II fit were set equal to the derived 6-values of Fe II in these components: 63 = 1.76 ± 0.14 km s~^ 
and 64 = 1.18 ± 0.11 km s"^ . 

JThe 6-values for component 2 in these ions were set equal to that derived for Fe II: b = 2.75 ± 0.07 km s~^. 
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Table 6. 


Component Fittin 


g Velocities 






Comp 2 


Comp 3 


Comp 4 


Comp 5 


Species 










X 


[km s-i] 


[km s-i] 


[km s~-^] 


[km s-i] 


CII 






41.25 ±0.12 


-28.0 ±0.05 


cir 


20.18 ±0.15 




41.7 ±0.3 




NI 


18.2 ±0.5 




41.10 ±0.05 




01 


20.38 ±0.12 


30.5 ±0.3 


41.20 ±0.03 




Mgl 


20.46 ± 0.05 


30.8 ±0.5 


41.20 ±0.21 




Mg II 


18.2 ± 1.0 


31.34 ±0.10 


41.20 ±0.03 


-28.37 ±0.19 


Al II 


19.88 ±0.23 


30.4 ±0.4 


40.9 ±0.3 




Si II 


20.1 ±0.3 


31.03 ±0.08 


41.20 ±0.05 


-30.2 ±0.8 


Si III 


21.6 ±0.5 


30.9 ±0.6 


41.1 ±2.0 


-28.3 ±0.4 


PII 


19.0 ± 1.2 








SII 


20.1*^ 








Cr II 


20.4 ±0.3 








Mn II 


20.4 ±0.3 


31.2 ±0.6 


40.2 ±0.8 




Fell 


21.16 ±0.03 


31.24 ±0.05 


41.21 ±0.05 


-30.1 ±0.9 


Ni II 


20.7 ±0.3 








Zn II 


20.6 ±0.6 









'^Central velocities in the LSR restframe and ±1(T errors (in km s~^) for the 
absorbing components 2-5. These values were derived through component fitting 
techniques as described in §3.3.2. 

'^The velocity for component 2 in the S II fit was held equal to that of Si II 
(see §3.3.2). 
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Table 7. Adopted Column Densities 









logA^x 


[cm -2] b 




Species 


logTVx" 


Comp 1 


Comp 2 


Comp 3 


Comp 4 




[cm -2] 


(3 km s-i) 


(20 km s-i) 


(31 km s-i) 


(41 km s-i) 


HI 


19.86 ±0.015^ 










CI 


12.96 ±0.03 


12.87 ±0.03 


12.21 ±0.09 






cr 


12.501°:!^ 


12.50t°;l° 








cr* 


< 12.16 


< 12.16 








CII 










13.75 ±0.05 


cir 






13.23 ±0.03 




12.36 ±0.05 


NI 






14.25 ± 0.08 




13.17±0.02 


I 


< 17.0'' 






13.58 ±0.03 


13.83 ±0.02 


Mgl 


12.61 ±0.02 


12.55 ±0.02 


11.72 ±0.02 


10.28 ±0.11 


10.52 ±0.04 


Mg II 


15.14 ±0.02 


15.08 ±0.02'' 


13.95 ±0.12 


12.50 ±0.02 


12.82 ±0.02 


Al II 






12.64 ±0.05 


11.47 ±0.05 


11.27 ±0.06 


Al III 


12.01 ±0.05 


12.01 ±0.05 








Si II 


15.14 ±0.02 


15.10 ±0.02'' 


13.98 ± 0.05 


12.78 ±0.02 


12.85 ±0.02 


Si ir 


11 47+0.12 


11 47+0.12 

-0.17 








Si III 






12.65 ±0.08 


11.83 ±0.10 


10.84t°:i? 


PII 


13.48 ± 0.02 


13.45 ± 0.02 


11.98 ±0.13 






SII 


15.21 ±0.02 


15.19 ± 0.02S 


13.82 ± 0.02 


< 12.85 


< 12.85 


S III 


13.82 ±0.02 


13.82 ±0.02 


< 12.59 


< 12.53 




Ca II'' 


12.35 ±0.06 


12.19 ±0.08 


11.85 ±0.08 




< 11.0 


Ti IP 


11.84 ±0.08 


11.78 ±0.09 


10.95 ±0.09 






Cr II 


12.62 ±0.03 


12.47 ±0.04 


12.03 ±0.05 






Mn II 


12.55 ±0.02 


12.48 ±0.02 


11.61 ±0.05 


10.82 ±0.09 


10.55 ±0.14 


Fe IP 


14.31 ±0.01'' 


14.13 ±0.02'^ 


13.78 ±0.02 


12.63 ± 0.02 


12.61 ±0.02 


Felll 


13.37^°:;? 


13.371°:?? 








Nill 


12.98 ± 0.06'^ 


12.86 ± 0.08'' 


12.35 ±0.04 
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Table 7 — Continued 



logNx [cm ^ 



Species 




Comp 1 Comp 2 Comp 3 Comp 4 

(3 km s"^) (20 km s'^) (31 km s"^) (41 km s'^) 



Cu II 
Zn II 



< 11.5 

12.57 ±0.08'^ 12.55 ±0.08"^ 11.16 ±0.08 



''Total sightline column densities derived from integrations of Na{v) profiles, unless 
otherwise noted. 

''Adopted column densities for the components considered here with ±la errors, given 
in atoms cm"^. Column densities for components 2-5 are derived from our component 
fitting analysis. Values for component 1 were derived through integrations of Na{v) 
profiles, unless otherwise noted. 

*^Derived using the continuum reconstruction method (see Appendix A). 
'^A 2a upper limit derived from the A1355.6 line. 
''Based upon the Mg II AA1239, 1240 lines. 
^Based upon the Si II A1808 line. 
SBased upon the S II A1250 line. 

^The Ca II data given here are from the profile-fitting results of Shull et al. (1977; 

see their Table 1). 

'The Ti II data given here are from the profile-fitting results of Welsh et al. (1997; 
see their Table 3). They have estimated the errors in the column density values to be 
~ 10% — 20%. The errors we have adopted in this table are slightly higher than this 
estimate. 

jWe have excluded the Fe II A2344 line and those lines with A < 1200 A in deriving 
this value due to possible oscillator strength uncertainties. 

'^The adopted column density given is derived from our component fitting 
measurements (see text). 



Table 8. CLOUDY Models of the fj, Columbae H II Region 



nil'' 
[cm-3] 




Rs' 
[pc] 








Log N{X^) [ 


cm ^J'^ 










S+ 


g+2 


P+ 


p+2 


N+ 


Si+ 




Mg+ 


Mn+ 


Fe+ 


0.02 


0.7 


275 


14.21 


13.79 


12.52 


12.06 


14.96 


14.53 


11.90 


14.47 


12.39 


14.18 


0.05 


0.5 


167 


14.23 


13.84 


12.53 


12.14 


14.99 


14.55 


11.95 


14.49 


12.48 


14.09 


0.2 


0.25 


84 


14.23 


13.84 


12.53 


12.14 


14.99 


14.55 


11.95 


14.49 


12.41 


14.20 


0.5 


0.15 


54 


14.22 


13.82 


12.52 


12.13 


14.98 


14.55 


11.80 


14.48 


12.40 


14.19 


01 


)S('r\('(l: 




ir).i<s 


13.82 


13.45 


< 12.(i 


15.7:'^ 


15.10 


12.1 5'' 


15.11 


12.48 


14.09 



''The assumed ambient hydrogen density (H atoms cm~^) for the H II region model. 
''The volume filling factor of the gas within the model H II region. 
'^The outer radius of the model H II region. 

'^Column densities of the indicated ions integrated through the model H II region, integrated from the 
star to Rs- Model calculations are not listed for Zn because of the large uncertainties in the atomic 
parameters for this element. 

"The N II result quoted here is from SY. The quoted uncertainties of this measurement are ±0.3, though 
given the large measured Wx, this value should be viewed as highly uncertain. 

^The Si IV result quoted here is from Brandt et al. (1998). The measurement is the weighted average of 
the lines AA 1394, 1403 A, log[A/'(Si+3)] = 12.16 ± 0.05. 
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Table 9. Normalized Gas-Phase Abundances 



Element 
X 


log(X/H)- 
+12.00 


[X/H]b 


w 


'/'SI'' 






Comp 1 
(3 km s-i) 


Comp 2 
(20 km s-i) 


Comp 3 
(31 km s-i) 


Comp 4 
(41 km s-i) 


C 


8.55 










-0.10 + 0.06 


N 


7.97 










-0.10 + 0.03 





8.87 








-0.52 + 0.04 


-0.34 + 0.03 


Mg 


7.58 


-0.30 ±0.02 


-0.42 + 0.03 


-0.18 + 0.12 


-0.31 + 0.03 


-0.06 + 0.03 


Al 


6.48 






-0.39 + 0.05 


-0.24 + 0.05 


-0.51 + 0.06 


Si 


7.55 


-0.26 ±0.02 


-0.37 + 0.03 


-0.12 + 0.06 






P 


5.57 


+0.05 ±0.02 


-0.04 + 0.03 


-0.14 + 0.13 






S 


7.27 


+0.08 + 0.02 










Ti 


4.93 


-0.95 + 0.09 


-1.07 + 0.05 


-0.53 + 0.09 






Cr 


5.68 


-0.92 + 0.03 


-1.13 + 0.04 


-0.20 + 0.05 






Mn 


5.53 


-0.84 + 0.02 


-0.97 + 0.03 


-0.47 + 0.05 


+0.06 + 0.09 


-0.28 + 0.14 


Fe 


7.51 


-1.06 + 0.02 


-1.30 + 0.03 


-0.28 + 0.02 


-0.11 + 0.03 


-0.20 + 0.03 


Ni 


6.25 


-1.13 + 0.06 


-1.31 + 0.08 


-0.45 + 0.04 






Zn 


4.65 


+0.06 + 0.08 


-0.02 + 0.08 


-0.04 + 0.08 







■^The logarithmic "solar" abundances of the elements, log(X/H)0, are used in deriving the gas normalized 

gas phase abundances. We have adopted the solar system mctcoritic abundances listed are from Anders & 
Grevessc (1989) except for C, N and O, which are photosphcric values from Grevesse & Noels (1993). 

''Wc present the sightline integrated values of [X/H] in this column. Thus [-^/H] = log{A^x/A^H} — 
log{X/H}0. Here we have used the value of logiVni = 19.86 + 0.015 (la systematic) derived in Appendix 
A as TVh- This treatment has neglected contributions from H+. However, our photoionization modelling 
(see §4) implies the corrections are relatively small (« —0.04 to —0.05 dex to the listed values). 

■^For the components 1 and 2 we have referenced the gas phase abundances to solar by comparing the 
column densities to that of sulfur. Thus [X/S] = \og{Nx/Ns} - log{X/S}0. 

'^For components 3 and 4 we have referenced the gas phase abundances to solar by comparing the column 
densities to that of silicon. Thus [X/Si] = log{iVx/iVsi} - log{X/Si}0. 
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Table 10. Dust-Phase Abundances for /j, Columbae Clouds 





(^/H)c^ 


106 (X/H)/ 


(X/Si)/ 


Element 




Comp. 1 


Comp 2 


Comp. 1 


Comp 2 




Solar Abundances (Anders & Grcvcssc 1989) 






Mg 


3.8 X 10-5 


23.2 ± 1.0 


13 ±8 


1.16 ±0.08 


1.5 ±1.2 


ol 


3.5 X 10-5 


20.0 ±1.0 


9±4 


1.00 


1.00 


Fe 


3.2 X 10-5 


30.80 ±0.08 


15.4 ± 1.1 


1.54 ±0.08 


1.8 ±0.8 


Ti 


8.5 X 10-8 


0.0777 ±0.0017 


0.060 ±0.006 


0.0039 ±0.0002 


0.007 ±0.003 


Cr 


4.8 X 10-^ 


0.442 ± 0.004 


0.18 ±0.04 


0.0221 ±0.0012 


0.021 ±0.011 


Mn 


3.4 X 10-^ 


0.302 ± 0.002 


0.224 ±0.015 


0.0151 ± 0.0008 


0.026 ±0.013 


Ni 


1.8 X 10-6 


1.690 ±0.019 


1.15 ±0.07 


0.085 ± 0.004 


0.13 ±0.07 




B-star Abundances (Gies & Lambert 1992)"^ 






Mg 


2.4 X 10-5 


12.0 ±0.8 


2±7 


0.58 ±0.05 


0.3 ±0.7 


Si 


3.3 X 10-5 


20.5 ±0.8 


10 ±3 


1.00 


1.00 


Fe 


5.2 X 10-5 


51.3 ±0.1 


38.0 ± 1.2 


2.5 ±0.10 


3.73 ±1.19 




B-star Abundances (Kilian-Montcnbruck et a 


■I. 19Mf 




Mg 


2.4 X 10^5 


16.6 ±0.5 


11±4 


2.05 ±0.15 


6±8 


Si 


1.6 X 10-5 


8.09 ±0.5 


1.7±2.0 


1.00 


1.00 


Fe 


2.5 X 10-5 


24.40 ± 0.004 


16.2 ±0.8 


3.01 ±0.19 


9± 11 



''Value for the "cosmic" abundance of each element using the reference abundances given in 
the table. 

''Dust-phase abundance of each listed element in parts per million. This value was derived 

by assuming S is undepleted and accurately traces the distribution of H I. In each case wc have 
used the (S/H)c value appropriate for the reference system. These arc (S/H)c = 1.9 x 10-5 
(Anders & Grevesse 1989); (S/H)^ = 1.6 x 10-5 (Gies & Lambert 1992); and (S/H)^ = 9.3 x 10-^ 
(Kihan-Montenbruck et al. 1994). 

'^Dust-phase abundance ratio of each listed element to Si, using the relative abundances given 
in the table. 

•^Gies & Lambert (1992) do not give abundances for Mg, Ti, Cr, Mn & Ni. We have adopted the 
Mg abundance derived by Kilian-Montenbruck et al. (1994) here and have not derived dust-phase 
abundances for these other elements. 

*^Kilian-Montenbruck et al. (1994) do not give abundances for Ti, Cr, Mn & Ni. We have not 
derived the dust-phase abundances for these elements. 
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Table 11. Mg, Fe, and Si Dust-Phase Abundances for Diffuse Clouds 



Cloud 106(Mg/H)d^ 106(Si/H)d* 106(Fe/H)d^ (Mg/Si)^'' (Fe/Si)/ [(Mg+Fe)/Si]/ Ref.'^ 



Warm Disk Clouds 



/U Col Comp. 1 


23.2±1.0 


20.0±1.0 


30.80±0.08 


1.16±0.08 


1.54±0.08 


2.70±0.12 


1 


HD 93521W.S " 


28±5 


22.3±1.8 


30.32±0.14 


1.24±0.23 


1.36±0.11 


2.6±0.3 


2 


C Oph Comp. A 


28.7±1.1 


24.5±1.0 


30.27±0.19 


1.17±0.07 


1.23±0.05 


2.40±0.08 


3 


HD 18100 


26.4±1.3 


20.0±1.7 


17.0±0.5 


1.32±0.13 


0.85±0.08 


2.2±0.5 


4 


HD 215733^ 


27.3±2.0 


24.8±1.5 


30.12± 0.21 


l.lOiO.lO 


1.22±0.07 


2.32±0.13 


5 


Weighted Mean: 


26.0±0.6 


22.4±0.6 


30.48±0.06 


1.17±0.04 


1.22±0.03 


2.46±0.06 




Warm Halo Clouds 


fj, Col Comp. 2 


13±8 


9±4 


15.4±1.1 


1.5±1.2 


1.8±0.8 


3.3±1.4 


1 


HD 93521w/'^ 


18±4 


15.5±1.8 


24.6±2.G 


1.1±0.3 


1.58±0.23 


2.7±0.4 


2 


HD 149881S 


22.5±1.4 


15.5±1.8 


26.5±0.5 


1.45±0.19 


1.70±0.20 


3.2±0.3 


6 


HD 215733 Comp. 6*^ 


26±7 


15±5 


27.23±1.06 


1.8±0.7 


1.9±0.6 


3.6±1.0 


5 


Weighted Mean: 


21.9±1.3 


15.3 ±1.2 


25.0±0.4 


1.37±0.16 


1.67±0.14 


3.11±0.20 





^Dust-phase abundance of the elements Mg, Si, and Fe in parts per million. These values were derived by assuming 
a species such as S or Zn is undepleted and accurately traces the distribution of H I. 

''Dust-phase abundance ratio of the elements Fe and Mg to Si. 
'^The sum of the dust-phase abundances of Fe and Mg relative to Si. 

•^Reference code for interstellar abundance data:(l) This work; (2) Spitzer & Fitzpatrick (1993); (3) Savage et al. 
(1992) and Sofia et al. (1994); (4) Savage k Sembach (1996); (5) Spitzer & Fitzpatrick (1995) (6) Fitzpatrick k Spitzer 
(1997). 

"^The data for HD 93521 were taken directly from Table 3 of Spitzer k Fitzpatrick (1993) and adjusted for differences 
in the Mg II oscillator strengths. The labels ws and wf refer to "warm slow" and "warm fast" clouds in their naming 
system. 

^Thc gas-phase abundance values used for deriving the dust-phase abundances for HD 215733 are based upon data from 
Table 4 of Fitzpatrick k Spitzer (1997). We have summed their components 4, 15, 17, and 19 for Mg, Si, Fe, and Zn. The 
dust-phase abundances above are based upon the resulting normalized gas-phase abundances [Mg/Zn]= —0.42 ± 0.08, 
[Si/Zn]= -0.39 ± 0.06, and [Fe/Zn]= -1.03 ± 0.04, using the suggested Mg II /-value of Fitzpatrick (1997). Following 

Fitzpatrick k Spitzer the dust-phase abundances quoted here were obtained by including a gas-phase abundance 
correction of —0.13 dex to account for a measured sub-solar abimdancc of Zn relative to S along this sightlino. 

SThe dust-phase abundances for HD 149881 are based upon data from Spitzer k Fitzpatrick (1995). These authors 
derive the gas-phase abundance of Mg for the integrated sightline ([Mg/Zn]= —0.39 ± 0.04 after scaling for the oscillator 
strengths used here). The Si and Fe values represent the sum of the column densities for their components 1, 2, 4, 6, 
and 8-10 and are normalized to Zn II in the same components using data from their Table 4. 

'^The data for this component of the HD 215733 come directly from Table 7 of Fitzpatrick k Spitzer (1997), with an 
adjustment for the [Mg/S] value to account for our adopted oscillator strength. 
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Table 12. GHRS Echelle Line Spread Functions at A1900 A 



Source 


Apert. 


Core 


Halo 




FWHM^ Frac.'^ 


FWHM*^ 


Frac.'' 


Robinson et al. 


LSA 


1.09 72.4% 


4.21 


27.6% 


This work 


LSA 


1.09 84.9% 


4.21 


15.1% 


Fitzpatrick & Spitzer'^ 


SSA 


0.95 95.9% 


5.00 


4.1% 



"^The FWHM in diodes of each component Gaussian. 

''The fraction of the total function in each component of the spread function. 

'^For the Fitzpatrick & Spitzer results we have exluded the very broad power- 
law wings in this treatment. Thus the fractions presented here are compared 
with the "core" plus near "halo" components. 
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Fig. 1. — Continuum- normalized intensity vs. LSR velocity absorption profiles for the species S II, P II, and 
Zn II. These species generally exhibit very low levels of depletion in the WNM. The stellar continuum for 
each line was approximated by a low-order polynomial (< 5) fitted to regions on either side of the line. For 
H Col, ulsr = fheiio — 19.9 km s^^. A telluric absorption line of O I in the P II A1301 profile is marked 
with the ® symbol. The S II A1253 observation contains a fixed pattern noise feature. We have marked 
observations made with the SSA. The absorption features extending from ~ 40 to 60 km s"-^ in the Zn II 
A2026 and P II A1301 panels are Mg I A2026 and O I A1302, respectively 
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Fig. 2. — As Figure but for the relatively lightly depleted species C II, O I, and N I. Telluric absorption 
lines of O I and N I are marked with the © symbol. 
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Fig. 3. — As Figure ^ but for the moderately depleted species Mg II and Si II. 
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Fig. 4. — As Figure P, but for the moderately to heavily depleted species Mn II, Cr II, Cu II, Al II, and Ni II. 
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Fig. 5. — As Figure nl but for the heavily depleted species Fe II. 
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Fig. 6.— As Figure 0, but for the non-dominant species Mg I, C IV, Si II, Si III, Si IV, Al III, S III, N II and 
Fe III. The C IV and Si IV profiles are from Brandt et al. (1998). The Fe III observations were made with 
the G140M grating. The N II are archival Copernicus data and were analyzed by ShuU & York (1977). 
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Fig. 7. — As Figure |l|, but showing transitions from excited states of Si II* A1264, C II* A1335, and several 
transitions of C I, CI*, and C I**. The term symbols for the lower level of the transition are given for each 
species. 
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Fig. 8. — Ground-based data for the /x Col sightline. The top panel is the WHAM spectrum of Ha towards 
Col taken with a 1° beam and 10 km s~^ resolution. An atmospheric OH line which coincides with 

component 4 is present at wlsr ~ +40 km s~^. The next panel shows the H I profile of Lockman, Hobbs, 
and ShuU (1986), taken with a 21' beam and 2 km resolution. The next three panels from top to bottom 
are the continuum-normalized profiles of Ca II and Na I from Hobbs (1978) and Ti II from Welsh et al. 
(1997). These data have velocity resolutions (FWHM) of Av w 4.5, 1.0, and 4.5 km s"""^, respectively. The 
data for Na I include contamination from telluric absorption lines. Those most likely to be present in the 
region of interest for this paper are near vlsr ~ —8, 5, and 14 km s~^ (see Hobbs 1978). 
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Fig. 9. Representative plots of the Na{v) profiles for the ions Mn II, Fe II, S II, and Mg II. For most of 
the absorbing components, these profiles exhibit little in the way of unresolved saturated structure. The 
Na{v) profiles of Fe II are discrepant near ulsr ~ +20 km s~^ (component 2), suggesting the presence of 
unresolved saturation in these transitions. 
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Fig. 10. — (a) Plot of the gas-phase abundances of various elements relative to S (see text). The x symbols 
represent the integrated abundances of the complex component 1 (t;LSR = +3 km s~^), while o symbols 
represent the adopted abundances for component 2 (ijlsr = +20 km s"-'^). (b) Same as (a) but with the 
observed gas-phase abundances of the Q Oph warm cloud included as the dotted line, and the spread of 
abundances for warm halo clouds (as taken from Sembach & Savage 1996) included as the hatched region. 
In both cases we have scaled the Mg II and Ni II abundances to be consistent with our adopted of the 
oscillator strengths (see text). 
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Fig. 11. — Normalized gas-phase abundances of several ionic species relative to S as a function of velocity. 
Plotted arc the ratios of P II, Zn II, Mg II, Si II, Mn II, Fc II, Ni II, and Cr II to S II, normahzcd to the 
solar system meteoritic abundances of Anders & Grevesse (1989), as a function of velocity. The error bars 
represent la uncertanties from photon statistics, continuum fits, and zero point uncertainties as well as a 
contribution for velocity shifts between the profiles. In the lower right of each panel is a bar representing 
the estimated size of the uncertainty due to contamination from material in an HII region surrounding fi 
Col. For Zn II and Mg II we have plotted both transitions since there are slight differences between them. 
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Fig. 12. — Plot of the predicted electron density, rie, as a function of temperature, Tg, using Mg I/Mg II and 
Ca Il/Ca III ionization equilibrium. The solid lines represent the trends in rig with Tg for the most likely 
ratios of these ionization stages. The dotted lines represent the la uncertainties (not including uncertainties 
in the atomic data). The values predicted by the intersection of the Mg I/Mg II and Ca Il/Ca III curves are 
(Te) = 6200 ± 300 K and (rie) = 0.27 ± 0.05 cm~3_ 
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Fig. 13. — Normalized gas-phase abundances of several ionic species in component 4 (t;LSR = +41 km s 
These species have been referenced to the column density of Si. 
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Fig. 14. — The ionic abundances of species relative to Sill and normalized by the X/Si ratio of the solar 
system as a function of ionization potential. The top plot shows the data for components 3 and 4 towards 
II Col. The bottom shows the same data with the data for 23 Ori (Trapero et al. 1996) and the VhSR = —58 
km s~^ component towards HD 72089 (Jenkins et al. 1998) overplotted as the solid and dot-dashed lines, 
respectively. The upper limits plotted for SlI are 2a limits. The ionization potential of Sill is 16.35 eV. 
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Fig. 15. — The square of the Doppcr parameter, 6^, derived for component 4 versus the inverse of the atomic 
mass, l/A, for the species N I, Mg II, Si II, and Fe II. The Hnc represents the best fit, yielding a temperature 
T = 4, 000 ± 700 K with a contribution from non-thermal velocities Vnt = 0.6 ± 0.4 km s~^. 
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Fig. 16. — Observed flux as a function of wavelength (LSR) in the G160M observations of the region 
containing interstellar Lya at A1215.67 A. These data have a signal to noise ratio of ~ 100 : 1. The tick 
marks the expected position of the line center. Absorption due to interstellar the Si III A1206.5 A transition 
is also visible in this spectrum. All of the broad features in this spectrum are due to stellar lines. 
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Fig. 17. — The best fit (filled circles) and ±2cr (points) values for the reconstructed continuum in the 
G160M and G140M observations (top and bottom, respectively). The best-fit column density of H I along 

this sightlinc is logA^(H I) = 19.87 ± 0.015. The thick line at the center of the line profile shows the range 
over which we have not reconstructed the profile, i.e., the range over which the signal to noise is less than 
10. 
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Fig. 18. — A comparison of the best-fit component models for the Mg II AA2796 and 2803 A transitions 
using the Robinson et al. LSF (1998; top row), that derived in this work (middle row), and the SSA data 
plus Spitzer & Fitzpatrick (1993) LSF. Also shown in the bottom row is the best fit model for Mg II A1239 A 
transition derived simultaneously with the stronger transitions using our derived LSF. The central velocities 
of the individual components are marked as ticks above the profiles. 
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Fig. 19. — A comparison of the instrumental LSFs for the GHRS at A1900 A. The top panel shows our 
derived LSF with the Robinson et al. (1998) LSF shown as the dash-dotted line, both for the LSA. The 
bottom panel shows our LSF compared with that derived for the SSA by Spitzer & Fitzpatrick (1993). The 
LSFs in a given plot are scaled so that each has equal area. 



